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A solution is obtained for the interior structure of an uncharged rotating black hole that accretes a 
collisionless fluid. The solution is conformally stationary, axisymmetric, and conformally separable, 
possessing a conformal Killing tensor. The solution holds approximately if the accretion rate is small 
but finite, becoming more accurate as the accretion rate tends to zero. Hyper-relativistic counter- 
streaming between collisionless ingoing and outgoing streams drives inflation at (just above) the 
inner horizon, followed by collapse. As ingoing and outgoing streams approach the inner horizon, 
they focus into twin narrow beams directed along the ingoing and outgoing principal null directions, 
regardless of the initial angular motions of the streams. The radial energy-momentum of the counter- 
streaming beams gravitationally accelerates the streams even faster along the principal directions, 
leading to exponential growth in the streaming density and pressure, and in the Weyl curvature 
and mass function. At exponentially large density and curvature, inflation stalls, and the spacetime 
collapses. As the spacetime collapses, the angular motions of the freely-falling streams grow. When 
the angular motion has become comparable to the radial motion, which happens when the conformal 
factor has shrunk to an exponentially tiny scale, conformal separability breaks down, and the solution 
fails. The condition of conformal separability prescribes the form of the ingoing and outgoing 
(f) ■ accretion flows incident on the inner horizon. The dominant radial part of the solution holds 

provided that the densities of ingoing and outgoing streams incident on the inner horizon are uniform, 
independent of latitude; that is, the accretion flow is "monopole." The sub-dominant angular part 
of the solution requires a special non-radial pattern of angular motion of streams incident on the 
inner horizon. The prescribed angular pattern cannot be achieved if the collisionless streams fall 
freely from outside the horizon, so the streams must be considered as delivered ad hoc to just above 
5-H , the inner horizon. 
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I. INTRODUCTION 



"No satisfactory interior solutions are known" §20.5] for rotating black holes. The purpose of this paper is to 
present nonlinear, dynamical solutions for the interior structure of a rotating black hole in the special case where it 
slowly accretes a collisionless fluid in a "conformally separable" fashion. A companion paper hereafter Paper 3, 
extends the solutions to charged rotating black holes. A concise derivation of the principal results of the present paper 
are given in in Paper 1 [|. A Mathematica notebook containing many details of the calculations is at [J] . 

Unlike the Schwarzschild @ geometry, the Kerr Q geometry contains an inner horizon as well as an outer 
horizon. In the classic analytically extended Kerr solution, the inner horizon is a gateway to delightful but unrealistic 
pathologies, including wormholes, white holes, timelike singularities, and closed timelike curves [8j. 

Sadly, the Kerr geometry fails at the inner horizon, because it is subject to the mass inflation instability discovered 
by Poisson & Israel (1990) 0. The inflationary instability is the nonlinear realization of the infinite blueshift at the 
inner horizon first pointed out by Penrose (1968) [lb1 |. Most studies of the inflationary instability have focussed on 
spherical charged black holes (see [ll| for a review). The first investigation of inflation inside rotating black holes 
was that of Barrabes, Israel & Poisson [l2j], who showed that when two light sheets pass through each other, a mass 
parameter defined by the product of the expansions along the two light bundles inflates. Ori [ill [TjJ and Brady and 
collaborators [HI, [l6| (see also references in [It}) considered inflation driven by a Price tail of ingoing and outgoing 
gravitational waves in the late time collapse of a rotating black hole. [l8[ found an exact mass inflation solution for 
a rotating black hole in 1+2 dimensions. 

Mass inflation is a generic classical mechanism that requires one essential ingredient to precipitate it: a source of 
ingoing (positive ene rgy) and outgoing (negative energy) streams near the inner horizon that can stream relativistically 
through each other [11| (in this context, positive and negative energy refer to minus the sign of the covariant t- 
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component pt of the momentum of a particle in a tetrad frame aligned with the principal frame, not to the conserved 
energy associated with translation invariance of the coordinate time t). Even the tiniest sources of ingoing and 
outgoing streams suffice to trigger inflation. As shown by [111, El, in spherical charged black holes, the smaller the 
streams, the more rapidly inflation exponentiates. This is a fierce and uncommon kind of instability, where the smaller 
the trigger, the more violent the reaction. The sensitivity of inflation to the smallest influence suggests that it would 
be difficult to avoid in a real black hole. 

Inflation is not particular about the origin of the ingoing and outgoing streams needed to trigger it: anything will 
do. Most of the literature has considered the situation where inflation is ignited by a Price [2CJ, |2l| tail of radiation 
generated during the initial collapse of the black hole. In real astronomical black holes, however, ongoing accretion 
of baryons and dark matter from outside probably soon overwhelms the initial Price tails. To illustrate how easy it 
is for accretion to generate both ingoing and outgoing streams at the inner horizon, consider a massive neutral dark 
matter particle that is slowly moving at infinity (energy per unit mass E — 1) in the equatorial plane of a Kerr black 
hole of mass M, and angular momentum parameter a. The dark matter particle can free-fall into the black hole if its 
specific angular momentum L lies in the interval (the following comes from solutions of the Hamilton- Jacobi equation 
derived in §IV Ajl 
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. (l + Vl + a/M.) <L< 2M. (l + y/1 - a/M.) . (1) 



The dark matter particle will become ingoing or outgoing at the inner horizon depending on whether its specific 
angular momentum L is lesser or greater than Lq given by 

U= l + j-ai/M? ' (2) 

which lies within the range (JT]) (hitting the upper limit of (fTJ in the case of an extremal black hole, a = M t ). In fact 
distant, slowly moving dark matter particles can become either ingoing or outgoing at the inner horizon provided that 
the inclination of their orbit to the equatorial plane is less than a value that varies from 90° (i.e. all inclinations are 
allowed) at a = 0, to 62° at a = 0.96M., to sin -1 ^/l/3 « 35° at the extremal limit a = M,. Thus cold dark matter 
particles falling from afar provide a natural continuing source of both ingoing and outgoing collisionless matter near 
the inner horizon. 

In black holes that continue to accrete, the generic outcome following inflation in spherical charged black holes is 
collapse to a spacelike singularity [ill HH [23| ■ 

Even if there were no accretion from outside, quantum mechanical pair creation would provide a source of ingoing 
and outgoing radiation near the inner horizon. Before the mechanism of mass inflation was discovered, [24| showed that 
electromagnetic pair creation in a spherical charged black hole would destroy the inner horizon. They suggested that 
pair creation, probably by electromagnetic rather than gravitational processes in view of the much greater strength 
of electromagnetism, would also destroy the inner horizon of a rotating black hole. [25[ considered a 2-dimensional 
dilaton model of gravity simple enough to allow a fully self-consistent treatment of pair creation and its nonlinear 
back reaction on the spacetime of a charged black hole. They found that pair creation led to mass inflation near the 
inner horizon, followed by collapse to a spacelike singularity. 

The inflationary instability means that the analytic extension of the Kerr geometry from the inner horizon inward 
never occurs in real black holes. Papers such as [26| that focus on the ring singularity of the Kerr geometry, while of 
interest in exploring how theories beyond general relativity might remove singularities, do not apply to real rotating 
black holes. 

Why does inflation take place at the inner but not outer horizon? The essential ingredient of inflation is the 
simultaneous presence of both ingoing and outgoing streams. Ingoing and outgoing streams tend to move towards the 
inner horizon, amplifying their counter-streaming. By contrast, outgoing streams tend to move away from the outer 
horizon, de-amplifying any counter-streaming. Classically, all matter at the outer horizon is ingoing: no outgoing 
matter can pass through the outer horizon. Outgoing modes will however be excited quantum mechanically near the 
outer horizon, leading to Hawking radiation. There have been various speculations that quantum effects could cause 
a quantum phase transition at the outer horizon [27l - l29l |. or prevent the outer horizon from forming in the first place 
[30(| . The present paper assumes that the outer horizon is not subject to a quantum instability. 

To arrive at a solution, this paper builds on physical insight gained from inflation in charged spherical black holes 
[Tlj [. Two fundamental features of inflation point the way forward. Firstly, inflation is ignited by hyper- relativistic 
counter-streaming between ingoing and outgoing beams just above the inner horizon. As shown in CVl regardless 
of their source or of their initial orbital parameters, near the inner horizon the counter-streaming beams become 
highly focussed along the ingoing and outgoing principal null directions. This has the consequence that the energy- 
momentum tensor of the beams takes a predictable form, Will The Kerr geometry is separable [31| , and the alignment 
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of the inflationary energy-momentum along the principal directions suggests that the geometry could continue to be 
separable during inflation. 

The second key fact is that inflation in spherical charged black holes produces a geometry that looks like a step- 
function, being close to the electrovac (Reissner-Nordstrom) solution above the inner horizon, then exponentiating 
to super-Planckian curvature and density over a tiny scale of length and time. This suggests generalizing the usual 
separability condition (|30j) to a more general condition (|46[) that departs, at least initially, by an amount that is tiny, 
but with finite derivatives. Unless one were specifically looking for step-like solutions, one would not think to consider 
such a generalization. 

The strategy of this paper is to seek the simplest inflationary solution for a rotating black hole. I assume that the 
black hole is accreting at a tiny (infinitesimal) rate, so that the Kerr solution applies accurately down to just above 
the inner horizon. This assumption is a good approximation for an astronomical black hole during most of its lifetime, 
since the accretion timescale is typically far longer than the characteristic light crossing time of a black hole. The 
approximation of small accretion rate breaks down during the initial collapse of the black hole from a stellar core, or 
during rare instances of high accretion, such as a black hole merger. 

For simplicity, I assume further that the geometry is axisymmetric, and possesses conformal time-translation 
symmetry (self-similarity), even though the accretion flow in real black holes is unlikely to be axisymmetric or 
self-similar. One might imagine that the assumption of conformal time-translation invariance would, in the limit of 
infinitesimal accretion rate, be equivalent to the assumption of stationarity, but this is false. As explained in §4.4 of 
(Tl| . and elaborated further in ^Illl of the present paper, the stationary approximation is equivalent to the assumption 
of symmetrically equal ingoing and outgoing streams at the inner horizon, whereas in reality the initial conditions of 
the accretion flow will generically lead to unequal streams near the inner horizon. The stationary approximation was 
first introduced to inflation by [32h34| , who called it the homogeneous approximation because the time direction t is 
spacelike inside the horizon. The present paper follow the convention of |35l . p. 203] in referring to time-translation 
symmetry as stationarity, even when the time direction is spacelike rather than timelike. 

I call the combination of conformal time-translation invariance coupled with the limit of small accretion rates 
"conformal stationarity," to distinguish it from strict stationarity. 

Motivated by the argument of the paragraph above beginning "To arrive at" , 1 pursue the hypothesis that the 
spacetime is conformally separable. By conformal separability is meant the conditions (|29|) on the parameters of 
the line-element §3§ and electromagnetic potential (|24[) that emerge from requiring that the equations of motion of 
massless particles be Hamilton- Jacobi separable, ^IV Al Whereas strict stationarity requires that the conformal factor 
take the separable form (|30p . conformal separability does not. Conformal separability implies the existence of a 
conformal Killing tensor, flVIl 

Finally, I assume that the streams that ignite and then drive inflation are freely-falling and collisionless. The 
present paper restricts to uncharged streams, while a companion paper i, Paper 3,' generalizes to charged streams. 
The assumption of collisionless flow is likely to break down when centre-of-mass energies between ingoing and outgoing 
particles exceed the Planck energy, but for simplicity I neglect any collisional processes. 



II. LINE ELEMENT 



Choose coordinates x^ = {x, t, y, </>} in which t is the conformal time with respect to which the spacetime is 
conformally time-translation symmetric (see i]IIII) . 4> is the azimuthal angle with respect to which the spacetime 
is axisymmetric, and x and y are radial and angular coordinates. Appendix [X] shows that under the conditions of 
conformal stationarity, axisymmetry, and conformal separability assumed in this paper, the line-element may be taken 
to be 



ds 2 



dx 2 



(dt — u! y d4>Y 



A„ 



(3) 



where 



v/1 



(4) 



The line-element is essentially the Hamilton- Jacobi and Schrodinger separable line-element given by [3ll . eq. (1)], 
except that the conformal factor p is left arbitrary, consistent with the weaker assumptions of conformal stationarity 
and conformal separability made here, as opposed to the stronger assumptions of strict stationarity and separability 
made by [31 1. 

Thanks to the invariant character of the coordinates t and <f>, the metric coefficients gu, gtrj>, and g^ all have a 
gauge-invariant significance. The determinant of the 2x2 submatrix of t—<p coefficents defines the radial and angular 
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horizon functions A x and A,, 



SttflVx/. - 9t<t> j A X A S 



(5) 



Horizons occur when one or other of the horizon functions A x and A y vanish. In physically relevant cases, horizons 
occur in the radial direction, where A^ vanishes. The focus of this paper is inflation, which occurs in a region just 
above the inner horizon, where the radial horizon function A x is negative and tending to zero. In this region the 
radial coordinate x is timelike, while the time coordinate t is spacelike. It is natural to choose the sign of the timelike 
radial coordinate x so that it increases inward, the direction of advancing proper time. 
Through the identity 



ds 2 



9fj,v dx^ 1 dx — ?7?rm^" fi 6 i/ dx^ dx , 



(6) 



the line-element ([3]) encodes not only a metric g^ Vl but a complete inverse vierbein e m fl , corresponding vierbein e m ' 1 , 
and orthonormal tetrad {7^, 74, -y y , 7^}, satisfying 7,„ -j n = r/ mn with rj mn the Minkowski metric. The tetrad defined 
by the line-element ([3]) aligns with the principal null tetrad, as will become evident from the fact that the Weyl tensor 
is diagonal in the tetrad frame (it has only spin-0 components, equation (|155[1 ). It is convenient to choose the time 
axis of the tetrad to lie in the radial x-direction, since that direction is timelike near the inner horizon. Explicitly, 
the inverse vierbein e m u is 



/ 1 



















V 



a 2 




(7) 



while the corresponding vierbein is 



m = _ 




(8) 



The convention in this paper and its companions is that dummy tetrad-frame indices are latin, while dummy 
coordinate- frame indices are greek. The tetrad-frame directed derivative is denoted d m , not to be confused with the 
coordinate- frame partial derivative d/dx^. Tetrad- frame and coordinate- frame derivatives are related by 



On 



dx^ 



(9) 



Whereas coordinate-frame derivatives d/dx^ commute, tetrad-frame derivatives d m do not. 



III. CONFORMAL STATIONARITY 



A fundamental simplifying assumption made in this paper is that the spacetime is conformally stationary, by which 
is meant that the spacetime is conformally time-translation invariant (that is, self-similar), and that the accretion rate 
is asymptotically tiny. Conformal stationarity generalizes the stationary (or homogeneous) approximation of [32l - [34j , 
which required equal ingoing and outgoing streams at the inner horizon, to the realistic case of unequal ingoing and 
outgoing streams. 

One might imagine that in the limit of asymptotically small accretion rate, the spacetime would automatically 
become stationary, d/dt = 0, but that is false. A characteristic of inflation is that the smaller the accretion rate, 
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the more rapidly inflation exponentiates . Even in the limit of infinitesimal accretion rate, the counter-streaming 
energy of ingoing and outgoing streams grows exponentially huge during inflation. The stationary assumption sets 
to zero some quantities that, although initially infinitesimal, nevertheless grow huge. As discussed in §4.4 of [ll[ and 
£|X Al of the present paper, the stationary assumption is equivalent to assuming equal ingoing and outgoing streams 
near the inner horizon. In reality, the relative fluxes of streams near the inner horizon depend on boundary conditions 
that generically do not lead to equal streams. 

In place of stationarity, a consistent approach is to impose conformal time-translation invariance, also known as 
self-similarity. Conformal time-translation symmetry allows the conformal factor p in the line-element @ to include 
a time-dependent factor, 

p = e vt p(x,y), (10) 

where the dimensionless factor p{x,y) is a function only of radius x and angle y (but not </>, given axisymmetry). 
Equation (I10[) says that if the conformal time t increases by an interval At, then the spacetime expands by a factor 
e vAt . The expansion is conformal, meaning that the spacetime keeps the same shape as it expands. The coefficient v 
can be thought of as a dimensionless measure of the rate at which the black hole is expanding. The conformal time 
t is dimensionless, and there is a gauge freedom in the choice of its scaling. A natural gauge choice is to match the 
change dt in the conformal time at some fixed conformal position well outside the outer horizon to the change d^KN 
in Kerr-Newman time measured in the natural units c = G = M» = 1 of the Kerr-Newman black hole: 

^KN ( } 

A/. V ' 

With that gauge choice, the accretion rate v is just equal to the dimensionless rate M, at which the mass M, of the 
black hole increases, as measured by a distant observer: 

dlnp d\nM. DM. ■ 

v = = — 7^r~ = t+ — = M » ■ 12 

at at otkn 

The accretion rate v is constant, so the mass of the black hole increases linearly with time as measured by a distant 
observer, M. = uiKN- 

The limit of small accretion rate is attained when 

v^O. (13) 

It might seem that the small accretion rate limit (fTBf would be equivalent to stationarity, but that is false. Whereas 
stationarity sets the accretion rate v to zero at the outset, conformal stationarity takes the limit (|13p after completing 
all requisite calculations, not before. 

In self-similar spacetimes, all quantites are proportional to some power of the time-dependent conformal factor e , 
and that power can be determined by dimensional analysis. The metric coefficients g^ v , inverse vierbein e m M , vierbein 
e m M , tetrad- frame connections Tfc/ m , tetrad-frame Riemann tensor Rkimm tetrad- frame electromagnetic potential A m , 
tetrad-frame electromagnetic field F mn , and tetrad- frame electromagnetic current j m , have respective dimensions 

SV ^ P 2 J e "V P J e ™ M K P 1 > ^klm OC p~ X , Rklmn OC p~ 2 , A m CX p° , F mn CX p" 1 , j m (X p~ 2 . (14) 

A quantity that is scale-free, such as the electromagnetic potential A m , is said to be dimensionless. 

The form of the time-dependent factor e vt in the conformal factor (fTQ| can be regarded as following from the fact 
that d m p must be dimensionless, which in turn requires that d\np/dt must be dimensionless, hence independent of t. 



IV. CONFORMAL SEPARABILITY 



Conformal separability is the proposition that the Hamilton- Jacobi equations of motion are separable for massless 
particles, but not necessarily for massive particles. Operationally, the proposition requires that the left hand side, but 
not the right hand side, of the Hamilton- Jacobi equation (|27|) be separable. A good part of mV Al and §JVB] below 
overlaps ground that is familiar since the work of [31| . These subsections are nevertheless needed to establish notation, 
to highlight where conformal stationarity and conformal separability differ from full stationarity and separability, and 
to provide the basis for the discussion in fjV] and for the derivation of the energy-momentum tensor of collisionless 
streams in ijVHI Subsection IIV CI covers new ground, showing that in the inflationary regime of interest here, even 
though the spacetime is only conformally separable, the equations of motion of massive particles are nevertheless 
Hamilton- Jacobi separable to an excellent approximation. Physically, massive particles are hyper-relativistic during 
inflation and collapse, and their trajectories are approximated accurately by those of massless particles. 
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A. Hamilton-Jacobi separation 



The equation of motion of a particle may be derived from Hamilton's equations, which express the derivative with 
respect to affine parameter A along the path of the particle of its coordinates x^ and associated generalized conjugate 
momenta tt^ in terms of its Hamiltonian H{x^, n^): 

~~dX ~ ~dx~» ' ~dX ~ ~£h~ ' ' ' 

The Hamiltonian of a test particle of rest mass m and charge q moving in a prescribed background with metric 
and electromagnetic potential is 

H = \g» V (jffi ~ qA^) {-K V - qA v ) . (16) 

The Hamilton-Jacobi method equates the Hamiltonian to minus the partial derivative of the action S with respect to 
affine parameter, H — —dS/dX, and replaces the generalized momenta with the partial derivatives of the action with 
respect to coordinates, 7r M = dS/dx 11 : 

i „„ fas . \ / os . \ os 



One integral of motion, associated with conservation of the rest mass m of the particle, follows from the fact that 
the Hamiltonian does not depend explicitly on the affine parameter. This implies that the Hamiltonian is itself a 
constant of motion: 

The normalization (|18[) of the Hamiltonian in terms of the rest mass m is equivalent to choosing the affine parameter 
A to be related to the proper time r along the path of the particle by 

dX = — . (19) 
m 

Two more integrals of motion follow from conformal time-translation symmetry, and axisymmetry. Conformal time- 
translation symmetry implies that the generalized momentum ~K t conjugate to conformal time t satisfies the equation 
of motion 

where the factor of — 2v comes from the time-dependent conformal factor e~ 2vt in the inverse metric g^ v in the 
Hamiltonian (|16[) . Equation (|2U|) integrates to 

m = — E — mvT . (21) 

In the small accretion rate limit v — > 0, this reduces to the usual conservation of energy, -K t = —E. The concern 
expressed in §1111 that some small quantities grow large during inflation does not apply here, because, as found in 
i jVIII Fl equation (|110[) . the proper time r experienced by a particle during inflation and collapse is always tiny, in 
the conformally stationary limit (this is checked explicitly at the end of mV C| . 

The two integrals of motion associated with conformal stationarity and axisymmetry correspond to conservation of 
energy E and azimuthal angular momentum L, 

dS dS T 

n = ^ = ~E, ^ = -=L. (22) 

Write the covariant tetrad- frame momentum pk of a particle in terms of a set of Hamilton-Jacobi parameters Pk , 

(23) 




Px Pt Py P<j> 



7 



and the covariant tetrad- frame electromagnetic potential Ak in terms of a set of Hamilton- Jacobi potentials Ak- 

Ax At Ay A(p I 




(24) 



u 



In Paper 3 it will be found that At is related to the enclosed electric charge within radius x, while A$ is related to 
the enclosed magnetic charge above latitude y. If magnetic charge does not exist, then A</> should vanish, but A,p is 
retained here to bring out the symmetry. The contravariant coordinate momenta dx K /dX = ek K p k are related to the 
Hamilton- Jacobi parameters Pk by 

dx K If P t UyP<j, u x P t P<l>\ 

-d\ = 7\ Px '^K' x + ~K y ~' Py '^K: + ~K y ')- (25) 

The tetrad- frame momenta pk are related to the generalized momenta tt k by pk = efc K 7r K — qAk, which implies that 
the Hamilton- Jacobi parameters Pk are related to the canonical momenta tt k by 

P x = -A x ir x - qA x , (26a) 
Pt = 7Tt + k^u x - qA t , (26b) 

Py = AylTy ~ qAy , (26c) 

P t jj = 7T0 + 7TtbJy - qAj, ■ (26d) 
In terms of the parameters the Hamilton- Jacobi equation (|17l) is 

p2 _ p2 p2 , p2 

y-+V = " mV ' (27) 

Separation of variables of the Hamilton- Jacobi equation proceeds by postulating that the action S separates as 
(equations (|28l) . (|29p . and below together constitute assumption III of [3l|, that the Hamilton- Jacobi equation 
separates "in the simplest possible way" ) 

S= ^m 2 \- Et + Lcf) + S x (x) +S y (y) , (28) 

where S x (x) and S y (y) are respectively functions only of x and y. The left hand side of the Hamilton- Jacobi 
equation (|2T|) is separable for arbitrary values of the constants E, L, and q provided that (cf. Appendix E|) 

uj x , A x , ,A X , At are functions of x only , , , 

D y , A y , Ay , A^ are functions of y only . ^ ' 

These are the conditions of conformal separability adopted in this paper. In the remainder of this paper, the black 
hole will be taken to be neutral, so that the electromagnetic potential A m is identically zero. The case of a charged 
black hole is addressed in Paper 3. 



B. Strict separability 

Full, or strict, separability, as opposed to just conformal separability, would require that not only the left hand side 
but also the right hand side of the Hamilton- Jacobi equation (|27p separates. This would require that the conformal 
factor p separates, as (this is eq. (43) of [3l|) 

P 2 = Ps =PI+P 2 y, (30) 

where p x is a function only of the radial coordinate x, and p y is a function only of the angular coordinate y. 
Equation ([3U)l holds during the electrovac phase prior to inflation, and it also holds during early inflation, when 
the conformal factor p remains at its electrovac value, but its radial derivatives dp/dx and d 2 p/dx 2 are becoming 
large. Equation (1301) breaks down as the conformal factor p begins to shrink from its electrovac value, presaging 
collapse. 

If p separates as equation (|3U)l . then the Hamilton- Jacobi equation (j2"T)) separates as 

p2 _ p2 \ p2 , p2 
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with K. a separation constant. The separated Hamilton- Jacobi equations (|31[) imply that 



P x = ±JP 2 - (AC + m>2) A, , (32a) 



P a = ±y/-P% + [K - m^pg) A, . (32b) 
The trajectory of a freely-falling particle follows from integrating dy/dx = —P y /P x , equivalent to the implicit equation 

dx dy 



p r ' (33) 

1 x 1 y 



The time and azimuthal coordinates t and <p along the trajectory are then obtained by quadratures: 

P f dx u>,,P<t, dy , Lj. T P t dx P^dy 
dt = -£-jr- + y * , dcj) = x 1 + -^-^ . (34) 

ix^x y y ixL^x lyL^y 



C. Conformal separability 

A feature of inflation, discussed in the next section, EfVl is that ingoing and outgoing streams of particles move 
hyper-relativistically relative to each other and to the no-going tetrad frame. The streams remain hyper-relativistic 
throughout inflation and subsequent collapse. One should not be too surprised that the trajectories of hyper-relativistic 
massive particles would be hardly distinguishable from those of massless particles. This subsection shows that 
Hamilton- Jacobi separability holds to an excellent approximation for massive as well as massless particles. The 
arguments are confirmed formally by showing that the difference (I39[) between the tetrad-frame momentum predicted 
by the Hamilton- Jacobi solution and the true momentum, integrated over the path of a particle during inflation and 
collapse, is adequately small. 

The Hamilton- Jacobi equation (|27|) can be separated as 

P x = Pt - [£ + m 2 (p 2 - p 2 )] A x , (35a) 
Py =-P^ + (IC- m 2 p 2 y ) A y , (35b) 

where K, is the same separation constant as before. The condition of separability is that the right hand side of 
equation (|35a[) is a function only of x, while the right hand side of equation (|35b|) is a function only of y. The latter 
follows from the conformal separability conditions (|29j) and the condition that p y is a function only of y, provided 
that TT t is treated as a constant. 

During the electrovac and early inflationary phases, p 2 equals its separable electrovac value p 2 , and equation (|35a[) 
simplifies to its electrovac form 

P 2 = P 2 -(IC + m 2 pl)A x , (36) 

whose right hand side is a function only of x, consistent with separability. During inflation and early collapse, the 
radial horizon function is tiny, |A X | < 1, so equation (|35a[) simplifies to 

Px = Pt , (37) 

whose right hand side is again a function only of x, consistent with separability. Once |A X | > f, equation (|57)l no 
longer holds, but during collapse the conformal factor p shrinks to a tiny value, with the net result that p 2 |A K | <C 1, 
so that equation (|35a[) simplifies to 

P 2 X = P 2 - {K- m 2 p 2 y )A x . (38) 

The term proportional to p 2 on the right hand side of equation (|3"5)) depends on y, apparently destroying separability. 
However, a feature of inflation and collapse, which will be discovered in jjVIII F[ is that the coordinates x and y of a 
freely-falling particle remain frozen at their inner horizon values throughout inflation and collapse. Thus p y remains 
constant along the trajectory of any particle, and the right hand side of equation (|38|) can be considered to be a 
function only of x, again consistent with separability of equation (|35al) . 

Equation (|38[) shows that during inflation and collapse the trajectory of a particle of rest mass m is accurately 
approximated by that of a massless particle with separation constant K-q = K. — m 2 p y . 



9 



That equations (|35|) provide accurate expressions for P x and P y for massive particles can be confirmed by considering 
the total derivative dpk/dx of the tetrad- frame momentum p^, equation (|23p . of a neutral particle of rest mass m along 
its trajectory. If the parameters P x and P y are taken to be given by equations (|35j) . with the parameters Pt and P^, 
from equations and 7r t constant (as opposed to from equation (l2~Tj) ). then the total derivative of the tetrad- frame 
momentum is 

dp k _ m 2 ( P y d{p 2 -pl) piyujyPA \ y/-A x y/A^ 1 m 2 i; f P t 1 P^uiy m y \ 

d\ ~ P 3 \2 dy + A y )\ P x ' ' P„ ' j + p \p,^A7 /=ST ^v 7 ^' V^J ' 

(39) 

The right hand side of equation (|39j) . which would be zero if the motion were exactly geodesic (and is in fact zero 
for massless particles, m — 0), is non-zero because the solution (|35p for P x and P y is not exact, for massive particles. 
The second of the two terms on the right hand side of equation (|39[) arises from approximating n t as a constant, and 
would disappear if irt were set to the more accurate value (|21[) . and r were replaced by its value as a function of 
respectively x and y in respectively Pt and P<j,. 

The deviation between the momentum p^ predicted by equations (|35p and the true momentum can be obtained by 
integrating equation (|39l) over the path of a particle during inflation and collapse. As shown in Appendix [D] the ratio 
Apk/pk of the deviation Apk = J (dpk / d\) dX to the momentum pk itself is of order ~ w 2 , equation (|D7|) . which may 
be considered adequately small. This confirms the accuracy of the Hamilton- Jacobi approximation (|35p . 



V. FOCUSSING ALONG THE PRINCIPAL DIRECTIONS 



A central feature of inflation, demonstrated in this section, is that as ingoing and outgoing streams approach the 
inner horizon, they see the opposite stream narrow into an increasingly intense, blueshifted beam focussed along the 
opposite principal null direction. The fact that near the inner horizon the energy-momenta of the streams becomes 
highly focussed along the principal null directions regardless of the initial conditions of the streams is what motivates 
the idea that the spacetime, which is separable in the Kerr-Newman geometry, may continue to be separable during 
inflation. 

A particle is said to be ingoing if P t , equation (|26bP . is negative, outgoing if P t is positive. Outside the outer 
horizon, P t is necessarily negative (ingoing) while P x can be either negative or positive. At the outer horizon, P t and 
P x are equal in magnitude, and continuous across the horizon. Inside the outer horizon, P t and P x switch roles: P x 
is necessarily negative (given that the sign of the timelike radial coordinate x is being chosen so that it increases as 
proper time advances, equation (|25p ). while P t can be either negative (ingoing) or positive (outgoing). A particle 
falling from outside the horizon necessarily has negative Pt as long as it is outside the horizon, but its Pt can change 
sign inside the horizon, if its angular momentum and/or charge are sufficiently large with the same sign as the black 
hole, as exampled in the Introduction in the paragraph containing equation (JTJ). 

A particle at rest in the tetrad frame has by definition tetrad-frame momentum p k = m{l, 0,0,0}, hence its 
Hamilton- Jacobi parameters Pk, equation (f23|) . are 



P k = -mp^A^l, 0,0,0} . (40) 

The tetrad rest frame defines a special frame, the no-going frame, where Pt = 0, at the boundary between ingoing 
and outgoing. 

The tetrad-frame 4-momentum p k of a freely falling particle, as seen in the no-going tetrad rest frame, is given by 
equation f|23[) . Near the inner horizon, where A x — > — 0, the no- going observer sees both ingoing and outgoing streams 
become hugely blueshifted and focussed along the ingoing and outgoing principal null directions: 

7£*- {l,±l,0,0}. (41) 

py/-A x 

Instead of the no-going observer, consider an ingoing or outgoing observer, of mass m', with Hamilton- Jacobi 
parameters Pj.. Near the inner horizon, where A x — > — 0, the ingoing or outgoing observer see particles in the opposite 
stream with hugely blueshifted 4-momentum p k (the following is equation (|4ip appropriately Lorentz-boosted in the 
radial direction), 

op' p 

P k ^ , * % r {l, ±1,0,0}, (42) 
m'p z (-A x ) 

in which the ± sign is + for ingoing observers and — for outgoing observers. An ingoing observer see a piercing beam 
of outgoing particles coming from the direction towards the black hole, focussed along the outgoing principal null 
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direction {1, 1, 0, 0}. An outgoing observer see a similarly intense beam of ingoing particles falling from the direction 
away from the black hole, focussed along the ingoing principal null direction {1, —1, 0, 0}. 

It is only particles on the opposing stream that appear highly beamed: particles in an observer's own stream appear 
normal, not beamed. The no-going observer is exceptional in being skewered from both directions, albeit with the 
square root of the energy and blueshift that an ingoing or outgoing observer experiences. 

VI. KILLING TENSOR 

Separability is associated with the existence of a Killing tensor, and conformal separability is associated with the 
existence of a conformal Killing tensor. 

A. Electrovac Killing tensor 



The separated Hamilton- Jacobi equation (|31[) can be written 

K mn PmPn = K , (43) 
where p m is the covariant tetrad-frame momentum (f!?3")) . and K mn is the tetrad- frame Killing tensor 

K mn = diagfog , -pi , p l , pi) . (44) 
The Killing tensor K mn satisfies Killing's equation 

D(kK r nn) = , (45) 

where denotes covariant differentiation, and parentheses denote symmetrization over enclosed indices. 

B. Early inflationary Killing tensor 

The separation of the conformal factor p as equation (|30p leads to the usual electrovac solutions [3l[ , but not to 
inflation. To admit inflation, it is necessary to go beyond equation (|30|) . This section proposes a modification of the 
conformal factor, and shows that there is an associated Killing tensor during early inflation. 

In spherically symmetric models, and in the limit of slow accretion, inflation has a step-function character at 
the inner horizon. This suggests generalizing the separation of the conformal factor p by allowing it to depart 
infinitcsimally from equation (|30[) . but with finite derivatives in the radial x direction. Specifically, the modified 
conformal factor is 

P = Pse vt ^ , (46) 

where p s is the usual separable conformal factor, equation (|30[) , and £, a function of x and y (not t or </>), is negligibly 
small, but with finite radial derivatives satisfying the hierarchy of inequalities 

d£ d 2 £ 

° ~ ^ ^ 'dx < da? ^ arly inflation ) ' ( 47 ) 



and with negligible angular derivatives, 



dy dy 2 ' ^ 



The inequalities (|47[) and (|48|l mean that £ is mainly a function of the radial coordinate x, its derivatives with respect 
to the angular coordinate y being small. The kind of function p that equation (|4"6"]l describes is one that takes a 
sharp turn, like a step function, in the radial direction. In Will CI it will be found that the initial conditions in the 
electrovac phase lead to a function £ that satisfies the conditions (|4"T)) and (l4"51) during early inflation. 

The separation (141)1) of the conformal factor p indeed proves to admit a Killing tensor, satisfying Killing's 
equation (|4"5)) . subject to the conditions (|47l) and (|48l) . The tetrad- frame Killing tensor K mn associated with the 
separation (|46l) is 

K mn = dmg(p 2 y , -pl , p 2 -p 2 y , p 2 - pl) , (49) 
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in which, to linear order in the small parameters v and £, 

p 2 -pl=pl + 2p 2 s (vt-0 ■ (50) 

Since v and £ are negligibly small, it might appear that the modified Killing tensor given by equation (|49[) is identical 
to the original tensor, equation (|44j) . Indeed, v can be set to zero in equation ()50[) without further delay. However, 
Killing's equation ([4"5"j) involves derivatives of the Killing tensor, which bring in non-negligible derivatives of £. Thus 
the modified Killing tensor (|49|) differs non-trivially from the original ([44]) . The Killing tensor (|49| satisfies Killing's 
equation (|45l) provided that not only t> and £ are negligible, but also d^/dy is negligible, but d^/dx may be large, 
consistent with conditions (|47|) and (j48|). 

The Killing tensor (l4l)l) applied to the tetrad-frame momentum given by equations ((25)) and ([55)1 gives 

K mn PmPn = K , (51) 

confirming that equations (1351) constitute a valid separation of the Hamilton- Jacobi equations for massive particles 
during early inflation. 



C. Conformal Killing tensor 

The previous two subsections, flVI Al and WIB1 showed that the spacetime possesses a Killing tensor in the electrovac 
and early inflationary stages, but not later. However, the spacetime possesses a conformal Killing tensor at all times, 
from electrovac through inflation and collapse. During early inflation, the inflationary exponent £ remains negligibly 
small while its radial derivatives grow large, conditions (|47l) . but during later inflation and collapse the inflationary 
exponent £ grows huge. 

The traceless part K mn of the Killing tensor (j49|) . 

K mn = K mn - \v" m K k k = |p 2 diag(l, -1, 1, 1) , (52) 

is a conformal Killing tensor [l], §35.3], satisfying 

D( k K mn ) - ^i](k m D l K n) i = . (53) 

The tensor K mn satisfies the condition (f53|) to be a conformal Killing tensor provided that the conformal separability 
conditions (|29j) hold, but without any restriction on the conformal factor p, and in particular without any restriction 
on the inflationary exponent £. 

In WIIIDI it will be found that during inflation and collapse the horizon function A x remains a function of radius 
x, as required by the conformal separability conditions f|29[) . only provided that the inflationary exponent £ is purely 
radial: 

£ is a function of x only . (54) 



VII. COLLISIONLESS FREELY-FALLING STREAMS 



The essential ingredient that triggers mass inflation is the presence near the inner horizon of ingoing and outgoing 
streams that can stream relativistically through each other. This paper adopts a general collisionless fluid as the 
source of energy-momentum that ignites and then drives inflation. In astronomical black holes, streams near the 
inner horizon will typically originate from accretion of baryons and cold dark matter. A combination of collisions and 
magnetohydrodynamic processes [III H3| are likely to keep baryons, electrons, and photons tightly coupled above the 
inner horizon, forcing them into a common ingoing or outgoing stream before inflation ignites. Dark matter, and also 
gravitational waves, which should behave like a collisionless fluid of gravitons in the high-frequency limit, can occupy 
the opposing stream, and stream relativistically through the baryonic stream without collisions, driving inflation. In 
the limit of small accretion rate considered here, the geometry above the inner horizon is accurately approximated 
by the electrovac (A-Kerr-Newman) solution, so the precise behaviour of the gas there is irrelevant. The strategy 
adopted in this paper and its companions is to seek solutions that hold from just above the inner horizon inward. 

The approximation of a collisionless fluid will probably break down when center-of-mass collision energies between 
ingoing and outgoing particles exceed the Planck energy. Such super-Planckian collisional processes, though probably 
important, are neglected in the present paper. 

It should be commented that general relativistic numerical treatments that model the energy-momentum as a single 
fluid with wave speed less than the speed of light are not satisfactory near the inner horizon, since such a fluid cannot 
support relativistic counter-streaming, and therefore artificially suppresses the inflation that would occur if even the 
tiniest admixture of an oppositely going fluid were admitted. 
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A. Occupation number 

The distribution of particles in a collisionless fluid is described by a scalar occupation number p) that specifics 
the number dN of particles at position x^ with tetrad-frame momentum p m = {p x ,p} in a Lorentz- invariant 6- 
dimensional volume of phase space, 

dN ='^$w- < 55 » 

Here d 3 x denotes the proper tetrad-frame 3-volume element measured by an observer at rest in the tetrad frame, not 
a coordinate-frame 3-volume element. The collisionless Boltzmann equation is 

where the affine parameter is dX = dr/m, with r the proper time of an observer at rest in the tetrad frame. The 
Boltzmann equation (151)1) asserts that the occupation number f(x^,p) is constant along phase-space trajectories. 
The tetrad-frame momentum 3-volume element d 3 p is related to the scalar 4- volume element d 4 p by 

( 2 " s » i ^' + '" 2 )W = *W' (57) 

where the Dirac delta-function enforces conservation of rest mass m. The Lorentz-invariant tetrad-frame momentum 
volume element from equation (|57[) translates into a 4-volume element of parameters P m with the Jacobian from the 
relation (|23[) between p m and P m : 



2 P *(2ttH) 3 v ' " \ p 2 A x p 2 A y J {2TTh)^p i A x A y 

This in turn translates into an element of orbital constants of motion, the energy E = —n t , angular momentum 
L = 7T0, and separation constant /C, with the Jacobian calculated from equations (|26bp . (|26d[) . amd ([35]): 

d 3 p a 2 dEdLdJC , , 

(59) 



2 P x (2ttH) 3 p 2 P x P y 4(27rft) s 



B. Number current 

The number density and flux of particles at any position form a tetrad- frame 4- vector n fc , 

Covariant number conservation follows from the collisionless Boltzmann equation l|56p . 

For a single stream with fixed constants of motion E, L, and /C, it follows from equation (|60|) and the constancy of 
the occupation number / along phase-space trajectories, equation (1561) . that the number current n k along the stream 
varies as 

2 

n k = Np k , N^-£—. (62) 

The denominators P x and P y would vanish where orbits turned around in radius x or angle y, and there would be 
cusps in the number current at such points. However, this never happens in the inflationary zone of interest in the 
present paper, because the radius x is a timelike coordinate, so P x never changes sign, while P y , given to an excellent 
approximation by equation (135b[) . remains essentially constant along any stream throughout inflation and collapse. 
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For massive particles, P x and P y along a stream are approximated accurately by equations (1351) . The accuracy of the 
approximations (I35[) can be checked by seeing how closely the covariant divergence D k n k that they predict vanishes. 
With P x and P y given by equations (|35|) . along with P t and P^ from equations (|26| and 7r t constant (as opposed 
satisfying equation (|21[) ). the covariant divergence of the single-stream number current n k given by equation (|62[) is 



D k n k 



m 2 N 




p 2 vP t 



(63) 



As expected, the divergence vanishes identically for massless particles, m = 0, but not for massive particles, because 
Hamilton- Jacobi separation is exact for massless particles, but not quite exact for massive particles. Since 



dlnN 1 fc „ fc 
~~dX~ = N DkU - DkP ' 



(64) 



and the momentum p k has already been checked, from equation (I39p . to be given accurately by expressions (|35[) for 
and P y , the deviation between the predicted and true logarithmic density In A can be estimated by integrating 
equation (|63[) over the path of a particle during inflation and collapse. As shown in Appendix [D] the resulting 
deviation A In A = J(d\nN/dX) dX is of order v 2 , equation (|D8|) . which may be considered adequately small. This 



again confirms the accuracy of the Hamilton- Jacobi approximation (|35[) . and the consequent expressions (1621) for the 
number density N and number current n k . 



C. Energy-momentum 



The tetrad-frame energy-momentum density T/-i of a system of freely-falling particles is 

d 3 p 
2pH2irh) 

Covariant energy-momentum conservation follows from the collisionless Boltzmann equation ([56 



T kl = I Pk pi f(x^p) ^7^3 . (65) 



D " T » = 1 ™ Dkf = l m fx = • (66) 

For a single stream with fixed constants of motion E, L, and IC, the energy-momentum is 

Tki = n k pi = Np k pi , (67) 

where the number current n k and number density N are given by equations (|62[) . Covariant energy-momentum 
conservation follows from number conservation and the geodesic equation, 

D k T kl = D k n kPl = piD k n k + = . (68) 

d\ 

In accordance with the expression (|23[) for the tetrad-frame momentum p k and the proportionality (|62p for N, the 
energy-momentum of a single stream is 

JfcJ = ; CX ; , (69) 

where 

A x for m = x , t , 

Aj, for m — y , (f> . ^ ' 

The behaviour of the collisionless energy-momentum will be explored in detail in §VHI[ but qualitative features of 
the behaviour are already apparent from equation (|69|) . As a stream approaches the inner horizon, A^ — > —0, the 
radial components of its energy-momentum grow large, because of the inverse factors of the radial horizon function 
in equation (|69[) . During inflation, the behaviour of the energy- momentum continues to be dominated by the radial 
horizon function A x , which is driven to an exponentially tiny value, causing the radial components of the energy- 
momentum to grow exponentially huge. During collapse, the conformal factor p shrinks, amplifying all components 
of the energy-momentum. 
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The trace of the energy-momentum tensor of a stream of particles of rest mass to, 

T% = Np k Pk = -Nm 2 , (71) 

can always be treated as negligibly small. During inflation, the trace is negligible because the incident accretion flow 
is negligibly small, in the conformally stationary limit. As inflation develops and collapse begins, the density N oc p~ 2 
increases as the conformal factor p shrinks, but the individual components of the energy-momentum increase more 
rapidly, as Tki oc p~ 4 , equation (|69|) . so the trace is never significant. 



VIII. INFLATIONARY SOLUTIONS 



This section presents the inflationary solutions that emerge from Einstein's equations with a collisionless source 
under the conditions of conformal stationarity, axisymmetry, and conformal separability assumed in this paper. 

Above the inner horizon, the A-Kerr-Newman electrovac geometry, Will Bl provides the background in which 
inflation ignites, Will CI As inflation develops, it back-reacts on the geometry, driving the conformal factor p and 
the radial horizon function A x from their electrovac forms. The equations governing the evolution of the conformal 
factor and radial and angular horizon functions are obtained in WIIIDI by separation of variables in two components 
of the Einstein tensor that have negligible collisionless source. The equations are solved to obtain the evolution of the 
conformal factor and radial horizon function in Will El The implications for the trajectories and densities of freely- 
falling streams are presented in §VIII Fl The solutions for the conformal factor and horizon functions are inserted 
into the remaining 8 Einstein components in Will Gl and in Will HI it is shown that the 8 Einstein components fit 
the form of the energy-momentum tensor of two collisionless streams, one ingoing and one outgoing. The solutions 
indicate that during collapse the angular (y and 4>) motion of the collisionless streams grows, and would begin to 
dominate once the radial horizon function is no longer small, \A X \ > 1. The last two subsections, Will II and Will .11 
address the effect of the angular motions to higher order, showing that, while the earlier results remain robust as long 
as angular motions are small, jA^j <C 1, the solutions fail when the angular motions become important, |A X | > 1. 



A. Charged black hole 



The case of a charged black hole is deferred to Paper 3, because the inclusion of electromagnetic currents, fields, and 
energy-momenta adds a whole extra layer of complexity to the solutions. Of course, a charged black hole is physically 
less interesting than the simpler case of an uncharged black hole, which is the focus of the present paper. 

Nevertheless, much of the results of this section and the next, Willi and EjlXI carry over essentially unchanged to 
the case of a charged black hole. In particular, the expressions (|88[) and (|124j) for the Einstein tensor are unchanged in 
the presence of an electromagnetic field produced by a collisionless source, and the solution, Will E[ for the evolution 
of the horizon function and conformal factor is unchanged, the only difference being that the derivative A' x of the 
electrovac radial horizon function at the inner horizon is altered by the presence of charge (because the electrovac 
horizon function (I77a[) depends on charge) , effectively changing the boundary conditions of the solution. The various 
differences between the charged and uncharged cases are detailed in Paper 3. 

To avoid repetition of the same results in Paper 3, the spacetime that provides the boundary conditions for the 
inflationary solution is referred to hereafter as "electrovac" rather than just "vacuum" . 



B. Electrovac initial conditions 



Electrovac solutions [l|,l31|, of which the physically relevant solutions are Kerr-Newman with a cosmological constant 
A, provide the boundary conditions for the inflationary solutions. 

The electrovac solutions are strictly stationary, satisfying v — 0, and strictly separable, so the conformal factor p is 
separable, equation (|30|) . Solution of the Einstein equations, Appendix [Cj leads to the standard results 



Ps = Pi + Pt = j^— c — \ tt , ; r , (72a) 




/ yu yx^x | yi yv^y r~7ou\ 

P* = \l 7i - i my? r~FTT\ ' Py = \l — nrz tti , . f ^ . ( 72b ) 
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and 

^7 = %V (To + /iwJ (.go - gi^x) , (73a) 
^ = 2^/(/i + /„*„) (, 9l - g Lj y ) , (73b) 

where / i A? Sen and <7i are constants set by boundary conditions. The sign of the square root for doj x /dx is the 
same as that for p x , while the sign of the square root for cku y /dy is the same as that for p y . Inflation leaves the 
separable factor p s in the conformal factor p, equation (|10[) . and the vierbein coefficients uj x and lo v unchanged from 
their electrovac values. 

For A-Kerr- Newman, the constants /o, /i, go, and gi are 

/o = , A = a' 1 ' 2 , g = a 3/2 , 9l = a 5/2 , (74) 

where a is the usual angular momentum parameter. The conformal factor and vierbein coefficients are given by 

p x — r — acot(ax) , p y — acosQ = —ay , (75) 



— . «,„ = asin 2 6> , a = 1 - uj x uj y = , R = V V 2 + a 2 . (76) 
K H 



The radial and angular horizon functions A x and A y are 

A 1 L 2M.r Q 2 m + Ql Ar 2 \ 

Ax B? y B? + B? 3 ) ' (77a) 

2 / Aa 2 cos 2 6»\ , 
A y = sm 2 eh + J , (77b) 

where M, is the black hole's mass, Q, and Q, are its electric and magnetic charge, and A is the cosmological constant. 
Inflation modifies the radial horizon function A x , but leaves the angular horizon function A y , along with uj x and ui y , 
unchanged. 

The full conformal factor p, equation (|4^]). involves an additional time-dependent factor of e vt (as well as an 
inflationary factor e~'> that equals unity away from the inner horizon). The parameters M t , Q,, Q,, a of the black 
hole coincide with the actual mass, charge, and specific angular momentum of the black hole at conformal time t = 0, 
and increase linearly with proper external time £kn- Physically, the cosmological constant A should not increase with 
time, but since it becomes completely overwhelmed during inflation by other exponentially growing energy-momenta, 
it can be included consistently in the description of the spacetime outside the inflationary regime (see the remarks in 
the paragraph following eq. 



C. Ignition 



This subsection outlines the behaviour of the inflationary exponent £ during the earliest phase of inflation, when 
the geometry is still electrovac, and collisionless streams are approaching the inner horizon. A more precise treatment 
that is valid throughout electrovac, inflation, and collapse starts in the next subsection, Will Dl 

As discussed in SjV] collisionless streams become highly focussed along the ingoing and outgoing principal null 
directions as they approach the inner horizon, causing the radial components T xx , T xtl and T tt of their energy- 
momentum to grow large. The Einstein combination G xx + G t t, which has zero electrovac source, is 



P 2 (G xx + Gtt) = - 2A 2 



dx 2 



c 2i 

dx 



- 2 



dH Ps /a) 

dx dx 



2v^ 
A x 



(78) 



The dominant term in this expression is the one proportional to the second derivative d 2 £/dx 2 of the inflationary 
exponent. Equating the Einstein component G xx + Gtt to the energy-momentum 8ir(T xx + T tt ) of collisionless streams, 
equation (|6"9")l . yields to leading order in 1/A.j. 



dx 2 



8tt£ NP 2 



(79) 
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where the sum is over collisionless streams. Because collisionless streams are hyper-relativistic near the inner horizon, 
the Hamilton- Jacobi parameters of every stream satisfy P t 2 = P x to an excellent approximation. The parameters N 
and P x of each component of the collisionless streams are sensibly constant as a function of radius x near the inner 
horizon preceeding inflation and during early inflation (but may be a function of angle y), and thus £ NP% is sensibly 
constant. The sum over streams £ NP% is of the order of the accretion rate v. 

In the situation of small accretion rate considered in this paper, the electrovac geometry provides an excellent 
approximation down to just above the inner horizon. Near the inner horizon, where A x — 0, the horizon function 
A x may be approximated by 



(x 



(80) 



where A' x = dA x /dx\ is the (positive) derivative of the electrovac horizon function at the electrovac inner horizon 
at x = X[ n . The derivative A^. is non-zero provided that the black hole is non-extremal, as should be true for any 
astronomically realistic black hole. Introduce the quantity u, a small positive parameter of order the accretion rate, 
u ~ v, defined by 



87r£iVPj 



(81) 



evaluated at the inner horizon x- ln . Later, A[ the overall accretion rates of ingoing and outgoing streams on to 
the inner horizon will be found to be proportional respectively to the combinations u =p v. Potentially u could be a 
function of angle y, but in WIIIDI it will be found that separability continues to hold as inflation develops only if u 
is independent of angle y. In terms of u, the second derivative ([75)1 of £ is 



9^ 

9a; 2 



u/A' x 



(*£ir. 



(82) 



Integrating equation (|82|) gives 



dx 



u/A' x 



u 

A7 



(83) 



where the constant of integration, established well outside the inner horizon, has been dropped because it is tiny 
compared to the retained term, which is diverging at the inner horizon x — > x m . Integrating equation (|83p in turn 
yields 



s U X in - 

t = -A> ln I "ST 



(84) 



where the constant of integration follows from requiring that £ starts at zero well outside the inner horizon, where 
x — > 0. Equation (|84[) shows that the inflationary exponent £ is the product of a small factor u ~ v and a term that 
diverges logarithmically at the inner horizon. Thus £ remains small even while its derivatives are becoming large. 
The inflationary exponent £ fulfills the conditions (j47|) and (|48|) anticipated in WIBI 



D. Equations governing evolution of the horizon function and conformal factor 

Inflation alters the spacetime geometry by changing the radial horizon function A^ and conformal factor p from 
their initial electrovac forms. Equations governing the evolution of the horizon function and conformal factor are 
obtained from the Einstein equations for the two diagonal components G xx — Gtt and G yy + G^. At least initially, 
these two components have negligible collisionless source in the conformally stationary limit. The angular component 
Gyy + G<f,$ has negligible collisionless source because the flow incident on the inner horizon has negligible energy- 
momentum, and inflation amplifies only radial components, not angular components. The component G xx — Gtt has 
negligible collisionless source because the trace of the energy-momentum of a collisionless source is always negligible, 
equation ([71]) . 

In Will Jl a non- vanishing collisionless source for G xx — G u and G yy + G^ will be taken into account, and it will 
be found that the results of the present subsection are robust. 
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Since the collisionless source is negligible for these two Einstein components, it is natural to seek homogeneous 
solutions of equations (|88p by separation. To achieve the desired separation, introduce U x and U y defined by 



U X = -^A X , 
Ox 



(85a) 
(85b) 



Initially, in the electrovac phase just above the inner horizon, U x equals the small parameter u defined by equation (|81l) , 

U x =u, (86) 

and U y is similarly small. As will be seen in Will El U x is driven by inflation to large values, but U y remains always 
small. Further, define X x , X y , Y x , and Y y by 



X x — 

Xy = 

Y x = 
Y» = 



dU* | 2 U *~ V2 
dx A x 



QUy 

dy A y 

dA T . d , 

— - + 3U X - A x — In 
ax ax 



it , t \ duJx 



(87a) 
(87b) 
(87c) 
(87d) 



In terms of U x , U y , X x , X y , Y x , and Y y , the Einstein components G xx — Gtt and G yy + are 



p (G xx ~~ Gu) — \ Ya 



P 2 ( G yy + G <t><t>) 



d In uj x y d In u y 
dx v dy 



-2X x + Y x —\n 
dx 



fo+flUx 



(gx-gou v ) duj y 



d In uj t d In lo„ 
Y — Y — 

dx y dy 



-2X y -Y y —\n 
dy 



dy 

h+fou v \ 



(g ~giuj x ) du x 



dx 



■*-y q — 

dy 



F^ln 
dy 



8Y X d 
X x + —^-Y x —ln 
ox dx 



duj y /dy 



U x (f + flLJ x ) 

duj x j dx 



Homogeneous solutions of these equations can be found by supposing that U x , X x , and Y x are all functions of radius 
x, while U y , X y , and Y y are all functions of angle y, and by separating each of the equations as 



1 ( foho+h 2 LJ x +fihiUjl fihi+h 2 ujy+foh ujy\ f h a +h 3 uj x fihi+hzU) y 



+ 



= 



(89) 



for some constants ho, hi, h%, and ^,3. If one attempts to separate equations (|88[) exactly, then the attempt fails 
unless U x and U y are identically zero, which is the usual electrovac case. But if U x is taken to be small but finite, 
then separation succeeds, and inflation emerges. If U x and U y on the second lines of equations (|88)) are treated as 
negligibly small, then separating the first lines of each of equations ([55]) according to the pattern of equation 
leads to the homogeneous solutions 



X x = 0, 

x y = o, 

(fa + fiu x )(ho + hiw x ) 



Y, 



Yy = 



dui x /dx 
(h + fo^ y )(hi + h UJ y ) 



duj y /dy 



(90a) 
(90b) 

(90c) 
(90d) 



If U x — U y — 0, then solution of the differential equations (I87c[) and (|87d|) with the homogenous solutions (|90c|) 
and (|90d[) for Y x and Y y , yields, subject to appropriate boundary conditions, the radial and angular horizon 
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functions A x and A y of the Kerr line- element. The separable solutions generalize to other electrovac spacetimes 
by admitting appropriate sources for Y x and Y y , Appendix [U] The solutions with a static radial electromagnetic 
field have a contribution G e mn = [{Q\ + Qt) / {p 2 pf)] diag(l, — 1, 1, 1), and those with a cosmological constant have 
Gmn = ~{p 2 /p 2 )^ J lmn- These electrovac contributions cease to describe a static radial electromagnetic field or 
cosmological constant when p ^ p s , but this happens only from the onset of collapse, by which time any electrovac 
contribution is overwhelmed by the collisionless energy- momentum, so the failure is unimportant (see §IV F of Paper 3 
for a more precise treatment). 

Solution of the angular behaviour during inflation is immediate. The vanishing, equation (|90b[) . of X y defined by 
equation (|87b[) implies that U y , which is initially neglible in the conformally stationary limit v — > 0, remains negligible 
throughout, and may be set to zero 

U y = . (91) 

The expression (|87d[) for Y y governing the angular horizon funtion A y is then unchanged from its electrovac form, 
and the angular horizon function A y thus retains its electrovac form during inflation and collapse. 

Down to just above the inner horizon, solution of the radial equations (|87a[) . (|87c[) . (|90a[) . and (I90c[) for X x and Y x 
leads to the usual electrovac form of the radial horizon function A^. But near the inner horizon, A x — > —0, the term 
proportional to \/A x on the right hand side of equation (|87a|) for X x starts to diverge, presaging inflation. In the 
vicinity of the inner horizon, where A x is near zero, expression (|87c[) for Y x simplifies to 

V, ^ • • (92) 
dx 

It will be found in the next subsection, §VIII El that the radius x remains frozen at its inner horizon value x- ln 
throughout inflation and subsequent collapse. Consequently u> x , hence Y x , are also frozen at their inner horizon values 
during inflation and collapse. Thus Y x is given by its electrovac value incident on the iner horizon, Y x = A' x where 
A' x = dA x /dx\ x , is the derivative of the electrovac horizon function at the inner horizon x = x- m . It follows that in 
the vicinity of the inner horizon the equations (|87a[) . (|87cp . (I90a[) . and (|90cp governing the evolution of U x and A^ 
reduce to 

, (93a) 
K ■ (93b) 

During the electrovac and earliest phase of inflation, when not only U x but also dU x /dx are negligible, which is true 
when \ A X \ > u, the Einstein components (|88[) separate without requiring that U x , which at this early stage equals it, 
be a function only of x. As inflation progresses however, continued separability, which requires that A^ be a function 
only of x, requires also that U x be a function only of x. It follows that the inflationary exponent £, given initially by 
equation (|84[) , must also be a function only of x. 

Below, equations f|l 24[) , it will be found that the energy-momenta of ingoing and outgoing collisionless streams are 
proportional respectively to U x =F v. If one or other stream vanished, then U x — v 2 would vanish, so equation (|93a[) 
would have no diverging term, and there would be no inflation. This is consistent with the physical argument that 
inflation requires the simultaneous presence of both ingoing and outgoing streams at the inner horizon. 

Equation (|93a[) indicates an instability only at the inner horizon, not the outer horizon. As streams approach the 
inner horizon, U x is driven away from zero because the horizon function is negative and tending to zero, A x — > — 0. 
By contrast, as infalling streams approach the outer horizon, the horizon function is positive and tending to zero, 
A x — > 0, causing U x to decay rather than grow. 



dU x 
dx 



dA x 

dx 



3C4 = 



E. Evolution of the horizon function and conformal factor during inflation and collapse 

This subsection integrates equations (|93[) . along with equation (|85al) . to determine the evolution of the radial horizon 
function A x and conformal factor p during inflation and collapse. 

The separation of the Einstein components ([55]) in the previous subsection, jjVIII D[ was premised on U x (and U y ) 
being negligibly small. However, the separation continues to remain valid during inflation and collapse when U x grows 
huge. The reason for this is that the dominant terms in the Einstein components (|88[) during inflation and collapse 
are of order U x /A x , coming from the expression (|87a[) for X x . Thus, once U x ceases to be negligible, the condition for 
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the validity of the separation becomes U x <C U 2 /\A X \, or equivalently |A X | <C U x . Thus the condition for the validity 
of the separation of the Einstein components (|88p is 

either U x < 1 or |A X | < U x . (94) 

It will be found in £ )VIII II that the conformally separable Einstein equations cease to be satisfied with a collisionless 
source once the angular motion of the collisionless streams becomes large, which happens when lA^I > 1. At this 
point U x is exponentially huge, equation (|107[) . Thus condition (|M)l remains well satisfied throughout inflation and 
collapse. 

The fact that the separation of the Einstein components (|88j) remains valid even when U x grows large, subject only 
to the condition is confirmed in WHIJi where the equations are solved to the next higher order in A x /U x . 
Integrating equation (|85a[) for U x and using equation (|93ap . gives the inflationary exponent £ as a function of U x : 

J A x J 2(U 2 ~v 2 ) A \u 2 ~v 2 ) ' y ' 

the constant of integration coming from £ = at U x = u, Equivalently, the inflationary part of the conformal 
factor p is 



\U 2 -v 2 



Inverting equation (|96|) gives U x in terms of £: 



U x = ^Jv 2 + (u 2 - v 2 )e 4 « . (97) 

Equation (I93b[) divided by equation (I93a[) yields an equation for the horizon function A x : 

dlnA x 3U X -A' 



dU x 2(U 2 - v 2 ) 



Equation (|98j) integrates to 

A x = - 



U 2 . 



3/4 



(U x + v)(u-v) 



{U x -v){u + v) 



K /(.*■») 



(98) 



(99) 



where the constant of integration is established by A^ ~ — 1 at U x = u. The precise constant is not important since 
near U x ~ u, the second factor on the right hand side of equation (|99[) is a number near unity taken to a large power 
A' x /(Av). The conformal factor p starts to depart significantly from its electrovac value p s when £ departs appreciably 
from zero, which occurs when U x is a factor somewhat greater than unity times u. At this point the horizon function 
is of order 

A, ~ e-V« , (100) 

which is exponentially tiny. 

The horizon function goes through an extremum, dA x /dU x — 0, where, according to equation 



which is of order unity. At horizon extremum, the inflationary exponent £ is 



2 

and consequently the inflationary part of the conformal factor p is 

/ v 1/2 

g _ ( 3yg - v 2 \ 1/2 



e = Ihif_^_) , ,i„2) 



e = A' ~ V ' (103) 



20 



which is becoming small. The horizon function at its extremum is 



3 / 2 /», „\ K/(^) 



eA' \ I u — v 



3Vu 2 - v 2 J \u + v 



e~ x l v , (104) 



which is exponentially tiny. 

The value of the radial coordinate x can be found by integrating equation (|93ap . 

x -^ = -J 2(U 2 -v 2 ) ■ (105) 

The integral can be expressed analytically as an incomplete beta-function, but the expression is not useful. Physically, 
equation (|105[) says that the radius x is frozen at its inner horizon value x- m during inflation and collapse, where U x 
is growing, while A x remains small. The radial coordinate x remains frozen even while the conformal factor p oc 
is shrinking. 

After the horizon function A x goes through its extremum, it starts increasing in absolute value as U x ^ 2 according 
to equation (IM|) . or equivalently as e 3 ^ according to equation (|97T) : 

A x « uV*Ao « e 3 « A , A = - (^) ^ , (106) 

where factors of order unity have been dropped compared to the exponentially huge factor Ao. In ijVIII II it will 
be found that the conformally separable Einstein equations can no longer be solved with a collisionless source once 
\^x\ ^ 1; the point at which the motion of freely- falling ingoing and outgoing streams become predominantly angular 
rather than radial. Equation (|106j) implies that at this point U x is 

CT, « |Ao 2/3 | atlA.lwl, (107) 

which is exponentially huge. In particular, the condition (|94p for the validity of the separation of Einstein components 
is well satisfied throughout inflation and collapse. 

During collapse, v and A^ in equations (|93j) may be neglected, and the equations yield a simplified equation for the 
evolution of the radius x, 

d\n(A x /U x ) _ U x ^ 
dx A x ' 

which integrates to 

x-x in = -^, (109) 

the constant of integration coming from x = x m at A^ = 0. Since \A X \/U X <C 1 throughout collapse, the radial 
coordinate x remains frozen at its inner horizon value x m to high precision through inflation and collapse. 



F. Trajectory and density of a freely-falling stream 

In the previous subsection, Will Ei it was found that inflation and collapse takes place over an extremely narrow 
zone in (conformal) radial coordinate x about the inner horizon value x- ln . Consequently the radial coordinate x 
attached to a freely-falling stream remains essentially frozen at its inner horizon value throughout inflation and 
collapse. The radial coordinate remains frozen even while the conformal factor p is collapsing to an exponentially tiny 
scale. 

The proper time r that elapses on the stream satisfies, equation (|25p . 



dr dX mp 2 (110) 

dx dx P T 



The right hand side of equation (|110[) is approximately constant during early inflation, while the conformal factor p 
is still close to its separable value p s , but then decreases as the conformal factor shrinks. Thus very little proper time 
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elapses on a stream during the entire of inflation and collapse. During collapse, the proper time r is even more frozen 
than the radial coordinate x, which itself is frozen. 

The angular coordinate y along the trajectory of the stream satisfies, equation (f33|) . 



dy 
dx 



P.. 



(Ill) 



The right hand side of equation (|111[) is a number of order unity as long as |A Z | <C 1. and then becomes less than 
unity when lA^I > 1. Thus the angular coordinate y is also frozen along the trajectory of a stream. It follows that 
quantities such as lo x and uj y , hence <r, are likewise frozen along the trajectory of a stream. 

The conformal time t and azimuthal angle (f> coordinates along the trajectory are given by equations (|34[) . which 
can be written 



dt 
dx 



1 

K 



Pi_ 

A* 



d(j) 

dx 



1 

K 



UyPt 



Pi, 



(112) 



As long as jA^I <C 1, the conformal time coordinate t along an ingoing (+) or outgoing (— ) stream is given by 



t = ± 



dx 



dU x 



= ± im 

4v 



(U x +v)(u-v) 



(U x -v)(u + v) 



(113) 



the constant of integration being established by t — at U x = u. 
factor p is 



(U x +v)(u- v) 



Consequently the time part e vt of the conformal 



(U x -v)(u + v) 



±1/4 



(114) 



Unlike x and y, the coordinate t is not frozen along the trajectory of the particle. Rather e vt , equation (|114[) . varies 
by a factor of order unity as U x increases from u to some large value. Once U x increases to some value much larger 
than u (such as U x ~ 1, since u <C 1), the time coordinate t freezes out. Once |A X | > 1, the relation t = ± f dx/ A x 
fails, but by that time t is already frozen out, so in practice equation (|114|) is valid accurately throughout inflation 
and collapse. 

The no-going tetrad frame has the special property that P t = 0. In the no-going frame the potentially large term 
proportional to Pt/A x in equation (|113|) vanishes, and t, like x, remains frozen throughout inflation and collapse. 
Physically, the factor e vt in equation (|114l) is the expansion factor of the black hole at the time that particles were 
accreted relative to the time that no-going particles were accreted. Equation (|114[) says that outgoing particles were 
accreted to the past of when no-going particles were accreted, when the black hole was smaller, while ingoing particles 
were accreted to the future of when no-going particles were accreted, when the black hole was larger. 

The furthest into the past that ingoing particles see (or into the future that outgoing particles see) is at the end 
of inflation when the geometry is collapsing. At this point the ratio of the sizes of the black hole at the times the 
ingoing (+) and outgoing (— ) particles were accreted is, from equation (|114[) with large U x , 



Paccrctc 



U + V 

u — V 



1/2 



(115) 



which is a number of order unity or a few. Equation (|115[) shows that what happens deep inside the black hole 
depends only on the finite past and future of the black hole, not on what happens at the initial moments of collapse, 
nor on the indefinite future. 

The density iV along a stream varies according to proportionality (|62j) . which given that a is frozen is 



TVcx 



1 



R 2 P X Py 



(116) 



The parameters P x and P y are given by equations (1351) . which are exact for massless particles, and whose accuracy 
for massive particles was established in ^IV CI 

As long as lA^I <^ 1, the Hamilton- Jacobi parameters P x and P y are constant along the trajectory of a stream, and 
the density iV ± of an ingoing (+) or outgoing (— ) stream simplifies to 



7V± oc _L oc e 2(?-0 ^ Ux T „ ) 



(117) 
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where the inflationary exponent £ and conformal time t have been eliminated in favour of U x using equations (|96[) and 
(|f f 4 [I . In this regime, the tetrad-frame momentum p^ is hyper- relativistic, and focussed along the radial direction, 

± l \ 1 1 1 P k / 11t n 

P * K pVtt' T tt'' iV '' i *J' " k= W^y- (U8) 

Note that pt, which are constant along the trajectory of the stream, are generically of order unity. 

It will be found in EIVIII Jl that the conformally separable Einstein equations cease to be satisfied by the energy- 
momentum of collisionless streams once the angular motions of the streams begin to exceed their radial motions, 
which happens when |A X | > 1. To see how this happens, it is necessary to (attempt to) follow the behaviour of 
collisionless streams into this regime. For lA^I > 1, as long as the spacetime remains conformally separable, the 
Hamilton- Jacobi parameters Pt, P y , and P^ all remain constant, but P x is no longer constant, varying in accordance 
with equation (I38|) . which can be written (the argument of the square root in the following equation is positive since 
A x is negative) 



Px 
P 

where p 2 + p 2 ^, a constant along the trajectory of the stream, is 



±y/l-(f4 + fil)A x , (119) 



„ „ K, — m 2 p 2 

t4+ >4 = — ■ ( 12 °) 



Putting together the dependence on p, equation (|117|) . and P x , equation (|119p . yields the density along an ingoing 
(+) or outgoing (— ) stream, 

N± ^ U xT v _ (121) 



i - Ok + n) A - 



The corresponding tetrad-frame momentum p^ is 



1 



whose angular components p y and p^ exceed the radial component p t once jA^j > 1. 



p " ocZ p\-- — 7=sT — ' t 7^a7' ^ h) ' (122) 



G. Einstein tensor 



The two Einstein components (|88|) led to equations for the evolution of the conformal factor and radial horizon 
function during inflation and collapse. The remaining 8 components of the tetrad-frame Einstein tensor Gki, which 
have zero electrovac source, may be written in terms of the expressions (|85|) for U x and U y , and (|87|) for X Xl X yi Y x 
and Y„: 



G, 



G t 
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„2 



i G x , 



p (G xy ± Gt y 
p 2 (G x0 ± G 
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" fcLx 7J 
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Y x ±v d 

3- ln 

— A T dx 
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(Uy =F IVUJy) 



(U x Tv) [A 
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- 2U, 
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(U x Tv) ( — %-^T2vu y ) 
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duj 



-A, 



dh\p 2 s VUJy dujy 

v dx a 2 dy 



„..-,, dlnp 2 
a dy ox 



X v T iv- 



U y du 

T 2 ' 

duj 



dy 



(123a) 
(123b) 

(123c) 
(123d) 



Most of the terms in equations (|123|) are negligible in the conformally stationary limit v — > 0. Terms are negligible 
because: (a) they are proportional to one of X x or X y , which vanish; (b) they are proportional to U y , which remains 
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negligibly small in the conformally stationary limit; or (c) they are proportional to v, and not proportional to inverse 
factors of A x , so remain negligible in the conformally stationary limit v — > 0. In addition, the term proportional to 
d\n(duj x /dx)/dx inside square brackets in equation (|123aj) may be neglected. This term might potentially become 
important when (A^l > 1, but it will be found below, equation (|143[) , that this term in any case disappears when the 
Einstein equations are solved to next order in A x /U x . 

With all negligible and sub-dominant terms discarded, equations (j 123[) simplify in the conformally stationary limit 

to 

p2 ( G** + Gu ± G ^ _ _L_ (JJ X qp „) (AL ± v) , (124a) 
p 2 (G xy ± Gty )=--y^=(U xT v)A y ^. , (124b) 
p 2 (G xct> ± Gn) = ± - \ = (U x T v) (^t*- T IVLJy) , (124c) 



y 

(?( Gm ~ G ** ±iG v ^=0, (124d) 

in which uj' x = duj x /dx\ x . and A' x = dA x /dx\ x . are the derivatives of uj x and the electrovac horizon function A x at 
the inner horizon x = x- m . The A' x in equation (|124aj) comes from Y x = A' x , equation (|93|) . The oj' x in equation ()124cj) 
comes from replacing dw x /dx, which remains frozen through inflation and collapse, by its inner horizon value u' x , as 
determined by equation (|73al) . 

Should not the factor A' x ±v on the right hand side of equation (|124aj) be replaced by A' x in the conformally stationary 
limit v — > 0, and likewise the factor A v Lo x /a 2 ^f2vujy on the right hand side of equation (|124cp be replaced by A y uj' x /a 2 l 
No. Each equation describes not one but two Einstein components, and the full expressions given are needed to capture 
both accurately. In the first case, the difference of equations (|124aP gives p 2 G x t = v(U x — A x )/(—A x ), and in the 
second case the sum of equations (|124c|) gives p 2 G xr / ) = —v(A y u)' x /cr 2 + 2ui y U x )/ W—A x A y . 



H. Streaming energy-momenta 

The form (|124j) of the Einstein components, derived under the conditions of conformal stationarity and conformal 
separability, fits to the form of the energy-momentum tensor T^i of a collisionless fluid consisting of two streams, one 
ingoing (+) and one outgoing (— ): 

Tkl=T+ + T-. (125) 
The energy-momenta of the ingoing (positive energy, p t < 0) and outgoing (negative energy, p t > 0) streams are 

Tti = ^Ptvt , (126) 



with densities 



and tetrad-frame momenta 



N ± = ^-(U xT v)(K±v) , (127) 



1 J 1 1 1 {A y d\np 2 /dy\ 1 / A y uj' x /(J 2 T 2ww : 



P 



V^A7 V^AI' A^± v J^^y A' x ±v 



Equations (|127l) and (| 128[) arc defined up to an arbitrary normalization that leaves T^i constant: the tetrad- frame 
momentum pk and density TV can be multiplied by some constant a and a~ 2 respectively. The corresponding Hamilton- 
Jacobi parameters , equation ([23| , of the collisionless streams are 



p±_J_i \d\np 2 Jdy _ A y u x /a 2 T2vu y \ 



(129) 



again up to an arbitrary normalization factor a. Equations (|127|) and (j 128[) agree with the behaviour (|117|) and (|118p 
of collisionless streams in the regime where the radial horizon function is small, 

\A X \ « 1 . (130) 
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Equations (| 127[) and (1128[) do not describe correctly the behaviour (] 1 2 1 [) and (1122[) of collisionless streams when 
the horizon function approaches and exceeds unity, lA^I > 1. Moreover the tetrad-frame momentum (|128p predicts 
that the purely angular components of the collisionless energy- momentum grow, becoming dominant when |A X | > 1, 
whereas it has been assumed from jjVIIIDI up to the present point that the angular components are negligible. The 
next two subsections, Will II and ^VIII J[ address these issues. 



I. Angular energy-momenta imposed by conformal separability 



As long as conformal separability holds, the angular components G yy — G^ and G yt p of the Einstein tensor must 
necessarily vanish, so that the 2x2 angular submatrix of the Einstein tensor must be isotropic (proportional to the 
2x2 unit matrix). The reason for this is that the expressions (|123d[) for those Einstein components are functions 
only of angle y, so the only kind of source of energy-momentum that they admit is one that depends only on y. By 
contrast, the densities N^, equation (|127|) . of the collisionless streams arising from inflation depend strongly on radius 
through U x . The only way that the collisionless streams can source the angular components of the Einstein tensor is 
that their angular energy-momentum tensor must be isotropic. 

Of course, as long as the angular components are sub-dominant, which is true when jA^I <C 1, there is no need for 
the angular components of the Einstein equations to be satisfied accurately, because the sub-dominant components 
have no effect on the remaining Einstein equations. This expectation is confirmed explicitly in the next subsection, 

gyiiui 

Nevertheless, as long as jA^I <C 1, an isotropic collisionless angular energy-momentum tensor can be accomplished, 
by allowing not just one ingoing and one outgoing stream, but several streams a, with densities 

Nt cx N* , (131) 

that sum to the totals prescribed by equation (|127[) 

J2K = n+, Y, N «= N ~- ( 132 ) 

a a 

Denote the means and mean squares of the angular components fi y and of the tetrad-frame momenta of the 
streams, equation (|118[) . by 

(^)^ EaN N i^ V , (M^ ^!^ , (133a) 

W = N+ + N- ' = N++N- ' (/W> = N++N- ' ( 33 } 

in which for the means the sum is over either the ingoing or outgoing stream, while for the mean squares the sum is 
over both ingoing and outgoing streams combined. The angular components \i y and /i^ must average to the values 
prescribed by equation (|128[) . 

(4) . f^) , (1 34a) 



The condition of angular isotropy requires that the mean squares of the angular components are equal, and the mean 
of their product is zero, 

</4> = (4) = ^ . (to**) = , (135) 
for some /x 2 . The angular components of the collisionless energy-momentum then form the isotropic matrix 



/PTyy = p'T^ = (N+ + N-)ff = UxA * V f . /r7-„„=l). 

on 



The Schwarz inequality, which states that the mean square of a distribution must exceed the squared mean, requires 
that 



2 ^ 
fi > max 



/ A+(/4)+A-( M -) \ 2 ( N+(»+)+N-(» 
V N+ + N- ) \ N+ + N- 



(137) 
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While contrived, there is no difficulty to construct a distribution of collisionless streams with the required mean 
momenta fll34[) and the conditions (|135l) and (|137[) on the mean squared momenta. 

On the other hand, the collisionless energy-momentum cannot be contrived to have isotropic angular components 
once |Aa;| > 1. In this regime, as long as conformal separability holds, the densities and tetrad- frame momenta 
of streams are given by equations (|121[) and (|122[) . The condition that the angular components be isotropic requires 
that, generalizing (|135|) . 

/ ^ =Y, N «-, ^ ( 138 ) 

If it happened that the means (|134j) of \i y and /x^, were the same, then the isotropy condition (|138[) could be satisfied 
for all values of the horizon function A x . But generically (e.g. for Kerr) the mean momenta differ, and the isotropy 
condition (I138[) on the squared momenta cannot hold for all A x . 



J. Inflation and collapse to next order 

In Will D[ equations governing the evolution of the horizon function and conformal factor were derived from the 
assumption that the Einstein components G xx — Gu and G yy + G<p<p had negligible collisionless source. It is true that 
the trace of the collisionless energy-momentum always remains negligible, equation (1711) . However, in Will HI it was 
found that the angular components of the collisionless energy-momentum, though initially negligible, grow, becoming 
dominant when (A^l > 1. In this subsection, the angular components are taken into account, which involves taking 
the Einstein equations to next order in A x /U x . It is found that the earlier results are robust as long as \A X \ <§; 1, but 
fail when angular motions become important, |A X | > 1. 

A collisionless source for the Einstein components G xx — Gu and G yy + G^ can be accommodated by admitting 
source terms for the quantities X x and Y x defined by equations (|87al) and (|87c|) . in much the same way that electrovac 
sources (a radial electromagnetic field and a cosmological constant) could be accommodated by admitting source 
terms for Y x and Y y . The sources for X x and Y x are conveniently written in terms of two arbitrary functions Fx and 
F Y : 

X 9 = -U m (p^l + F x ^ , (139a) 
Yx = A x ( g^ft + 2 (/°fli+Ago) _ Fy ) . (139b) 



dx du x /da 

During early inflation, when both U x and A^ are negligibly small, both sets of source terms are negligible. The source 
terms are of order A x /U x compared to the principal terms in X x and Y x , and remain sub-dominant through inflation 
and collapse. The contributions d\n(p s /a 2 )/dx and dlnp s /dx to the source terms depend on angle y as well as radius 
x, breaking separability, but the breakdown is unimportant as long as condition (|94j) holds, which is well satisfied 
through inflation and collapse. 

It is convenient to define new quantities X x and Y x that concatenate X x and Y x with their source terms: 

Y, s ^ + W, + A. (p^§p- + F r ). (140b) 

The equations governing the evolution of U x and A x are then 

X x = , (141a) 
Y X =A' X , (141b) 

generalizing equations (|93[) . Equations (|141l) constitute the evolution equations for U x and A^ taken to next order in 
A x /U x . As long the condition for the validity of the earlier evolution equations holds, the evolution of U x and A x 
is essentially unchanged, and all the results of Will El carry through unchanged. This is as expected physically: the 
evolution of the conformal factor and horizon function should be essentially unaffected by sub-dominant contributions 
to the energy-momentum. 
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Expressions for the Einstein components G xx — Gtt and G yy + in terms of the higher order X x and Y x defined 
by equations (|139al) and (|139h>[) are 



p (G xx — Gtt) — — %X X — 2Y X 



- din p., 



i).r ' ^ y dy ^ v dy^ U (er 4 dy ^ X 



2A 7 



1 / 1 doj, 



4 \a 2 dy 



F Y 



dlnp s 

dx 



dy 



a 20 \ dy 



(142a) 



P 2 (Gyy + G$$) = — 2X y + 2Y f 



d In p s 
dy 



X x 



dx 



Y x 



A 3 . , 

4 \a 2 dy 

while the Einstein component (|123al) becomes 



3/1 duj,, 



dF Y 
dx 



F Y + F Y 



dH Ps /a 2 ) 
dx 

dln(p 2 /a 2 
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dx 
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+ U x (- F x + 3F Y ) 
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(142b) 



(143) 



with expressions (I123bl) . (|123cp , and (I123dl) remaining unchanged. 

During collapse, when U x is large, the dominant terms contributing to the right hand sides of equations (|142j) are 
those proportional to U x . Expressions for the functions Fx and F Y follow from requiring that the [/^-dependent 
terms on the right hand sides of equations (1142)) are, in accordance with Einstein's equations, equal to 87r times the 
corresponding collisionless energy-momenta. For lA^I 1, the relevant collisionless energy-momenta are the isotropic 
angular components (|136[) . and a negligible collisionless trace, Tj? = 0. The resulting functions Fx and F Y are 



F 



x 



(144) 



The functions are of order unity, and thus, as anticipated, their inclusion has negligible effect on the evolution 
equations (|141[) as long as condition (|94|) is true, which as shown in Will El is well satisfied through inflation and 
collapse. 

If the functions Fx and F Y from equations (|144[) are inserted into the right hand side of equation (1143[) . then the 
result agrees with the earlier expression (|123ap except that the sub-dominant term proportional to dln(duj x /dx)/dx 
in the earlier expression disappears. The net result is that, when the conformally stationary limit is taken, the same 
set of expressions (|124[) is obtained for the Einstein components as found previously. This confirms that taking into 
account the purely angular components of the collisionless energy-momentum has essentially no effect on the solution, 
as long |A X | 1. 

On the other hand, as already discussed in EIVHI I| once lA^I > 1, generically (e.g. for Kerr) the angular components 
of the collisionless energy- momentum can no longer be arranged to be isotropic, equation (|138l) . so the conformally 
separable Einstein equations cannot be satisfied with a collisionless source. Even if it could be arranged that the 
angular components were isotropic, once lA^I > 1 the situation is further complicated by the fact that dio x /dx, which 
up to this point has been frozen at its electrovac value, must also change in order to satisfy the Einstein equation for 
Gtcj,- Further discussion is unwarranted in this paper. 

In conclusion, the solution found earlier in this section, up to ijVIII HI holds as long as |A X | <C 1, but breaks down 
when the angular motions of the collisionless streams exceed their radial motions, which happens when lA^I > 1. 
What happens after angular motions become important is undetermined. 



IX. MASS AND CURVATURE 



A. Mass inflation 



The term mass inflation comes from the fact that in charged spherically symmetric models of black holes, the 
interior mass, or Misner-Sharp (38| mass, inflates exponentially. In spherical black holes, the interior mass M can be 
defined by 



2M 
r 



1 



-8 m r d m r , 



(145) 
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where r is the circumferential radius, which plays the role of a conformal factor, and d m is the directed derivative in 
any tetrad frame. An analogous scalar quantity in the rotating black holes considered in the present paper is 



2M 
P 



l = M=M 2 



-d m pd mP 



where M x and M. y are dimensionless radial and angular mass parameters defined by 
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(147a) 
(147b) 



The mass parameter proposed by [12| for rotating black holes is the radial mass parameter M. x . In the Kerr-Newman 
geometry (A = 0), the interior mass M defined by equation () 146(1 goes over to the black hole mass M. far from the 
black hole, 

M -> M. as r -> oo . (148) 
In terms of the quantities U x and U y defined by equations (|85[) . the dimensionless mass parameters are 
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During the electrovac phase prior to inflation, the radial mass parameter M. x is proportional to the horizon function 
A x , 



M x = 



dm p s 

dx 



A,. 



(150) 



which approaches some small value near the inner horizon A x — > —0. The radial mass parameter M. x given by 
equation (|149a|) reaches a minimum, signalling the start of inflation, when (here U x equals its initial value u) 



d_\n£s 
dx 



A x = - 



(151) 



At this point the streaming energy and pressure in the no-going tetrad frame are comparable to unity in natural black 
hole units, c — G — M, = 1, while the mass parameter is small, of order M. x ~ v. Once inflation gets going, the mass 
parameter is 



U 2 



(152) 



which can be recognized as the principal term driving the evolution of U x , equation (I93ap . During inflation the radial 
mass parameter increases exponentially. 

During inflation and collapse, the ratio of the mass parameter M. x to the streaming energy T xx in the no-going 
frame is 



M x 
8irT xa 
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2 u x 



U X A' X 
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(153) 



which remains always of order v, increasing mildly as U x increases from u to some large value. 

3 /2 

After the horizon function has gone through its extremum, it increases in absolute value as A x oc U x , so the 

1 /2 

dimensionless mass parameter varies as M. x oc U x , or equivalcntly as M. x oc 1/p. The interior mass M defined by 
equation (|146p thus goes to a (huge) constant during collapse, consistent with behaviour of the interior mass during 
collapse in charged spherical black holes, §4.3 of fu\ . 
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B. Weyl curvature 



The only non-vanishing component of the Weyl tensor is the complex spin-0 component. The fact that only the 
spin-0 component is non-zero defines the spacetime as Petrov type D. Subject to the conditions of conformal time 
invariance (not necessarily conformal stationarity) and conformal separability assumed in this paper, the polar (real) 
part of spin-0 Weyl component, equations (|B5[) - (IB7[) . may be written in terms of the quantities U y , X y , and Y y defined 
by equations (|85b[) . (|87bl) . and (|87d[) . as well as the mass parameter M. defined by equation (I146p . as 
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which is valid throughout the electrovac, inflationary, and collapse regimes. Since U y and X y are negligible in the 
conformally stationary limit, the polar spin-0 Weyl component reduces to 
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(155) 



I have not been able to find any enlightening expression for the axial (imaginary) spin-0 component of Weyl tensor, 
beyond that already given as equation (|B8[) . 

During inflation and collapse, the Weyl curvature is dominated by the mass term, 



(156) 



X. BOUNDARY CONDITIONS 



In ^ Villi it was found that a conformally stationary, conformally separable solution exists for the interior structure 
of a rotating black hole, from electrovac through inflation to collapse, in which the Einstein equations are sourced 
by the energy-momentum of ingoing and outgoing collisionless streams. Because the accretion rate is asymptotically 
tiny, the streams have negligible energy-momentum above the inner horizon, and therefore have no effect on the 
geometry above the inner horizon. From the point of view of boundary conditions, what is important is the form of 
the collisionless streams incident on the inner horizon. 



A. Density of collisionless streams incident on the inner horizon 



The most important boundary conditions are those on the radial (x-t) components of the collisionless energy- 
momentum incident on the inner horizon. The radial components of the energy-momentum of ingoing (+) and 
outgoing (— ) collisionless during inflation and collapse are, equations (j 127[) and (|128|) . 

N ± 

P \&x\ 

The initial values of these components are set by the densities N^, equation (|127|) . of the ingoing and outgoing 
collisionless streams incident on the inner horizon, which are, since U x = u initially, 

N ± = -±-(u T v)(A' x ±v) . (158) 

The densities N ± of incident ingoing and outgoing streams are proportional to u =F v respectively. Both ingoing and 
outgoing streams must be present for inflation to occur, so both densities must be strictly positive. Moreover, at least 
classically, the black hole must expand as it accretes, so v must be positive. Thus the accretion rates u and v, both 
of which are to be considered small, must satisfy 



u > v > 



(159) 
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Positive v implies that in the no-going tetrad frame the density of the outgoing stream exceeds that of the ingoing 
stream during inflation and collapse. Equivalently, the center-of-mass frame is outgoing. An outgoing density 
that exceeds the ingoing density ensures that the black hole's angular momentum, which in self-similar solutions 
is determined by the angular momentum of the accreted streams, is positive. The case v = corresponds to equal 
ingoing and outoing streams, which is the stationary (or homogeneous) approximation first applied to inflation by 

Equation (|158[) prescribes that the densities of incident ingoing and outgoing streams must be uniform, 
independent of latitude y. In other words, the accretion flow must be "monopole" in density. It makes physical 
sense that the condition of conformal separability would require a high degree of symmetry of the incident accretion 
flow. 

The small accretion rate parameters u and v completely characterize the inflationary solution. The evolution of the 
inflationary exponent £ and horizon function A x found in Will El is determined entirely by these parameters, along 
with the derivative A' x of the electrovac horizon function at the inner horizon. 

The positivity of the ingoing and outgoing densities N ± , equation (|158|) , requires also that A' x ± v be positive. This 
means that the inflationary solutions do not apply to extremal black holes, whose inner horizons coincide with their 
outer horizons, and for which A' x = 0. 



B. Angular motion of collisionless streams incident on the inner horizon 

The angular components of the momenta of the collisionless streams are sub-dominant during inflation and collapse. 
As seen in EIVIII J[ the evolution of the inflationary exponent £ and horizon function A x , and of the radial components 
of the collisionless energy-momentum, are unaffected by angular motions until angular motions become important, at 
which point the solution fails, which happens when the geometry has collapsed to exponentially tiny scale. As far as the 
radial solution is concerned, a sufficient condition for conformal separability to hold is the condition (I158[) on the density 
of the accretion flow. If however the Einstein equations are required to hold also for the sub-dominant radial-angular 
components of the collisionless energy-momentum, which is a more stringent constraint on conformal separability, 
then the angular components of the tetrad-frame momenta of the ingoing and outgoing collisionless streams must be 
as given by equation (I128|) . The corresponding Hamilton- Jacobi parameters P^ are given by equation (I129p . The 
condition (|129[) prescribes that the angular components of the Hamilton- Jacobi parameters Pj^ do not vanish, but 
rather vary in the given fashion with latitude y. 

Equation (|129[) requires that the collisionless streams have some net motion in the angular y direction. One might 
perhaps have anticipated that conformally stationary flow might require that each particle in the collisionless stream 
fall along a surface of constant latitude y, which would happen if P y for each particle were identically zero. Trajectories 
at constant y occur when P y is not only zero but also an extremum with respect to variation of x or y at fixed constants 
of motion E, L, K.. In fact equation (|129j) shows that the required P y is not zero. 



C. The angular boundary conditions cannot be achieved with collisionless streams accreted from outside 

the outer horizon 

Can the relations (|129[) governing the angular motions of streams incident on the inner horizon be accomplished by 
collisionless streams that accrete from outside the horizon? As will now be shown, the answer is no. 

Equation (|129[) shows that the ratio Pf 1 /P^ incident on the inner horizon for each of the ingoing (+) and outgoing 
(— ) streams is 

If = A' x ±v 

If the ingoing and outgoing streams contain multiple components, then the ratio (|160[) is a density-weighted average. 
If one insists that every particle originate from outside the outer horizon, then P t must be negative at the outer 
horizon for every particle. The relations (|26[) with vanishing electromagnetic potential (the case of a charged black 
hole is deferred to Paper 3) then imply that Pt and of the particle at the inner horizon must satisfy the inequality 

P t < P$ ~ Ux ' out , (161) 

1 ^x^out^y ,in 

where the subscripts in and out denote values respectively at the inner and outer horizon. The condition (|161[) 
excludes half of the Pt-P<j> plane. As illustrated in Figure [TJ the right hand sides of equations (I160p and (|16ip have 
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separability condition 



-1.0 -.5 .0 .5 1.0 

cos 6) 

FIG. 1: Comparison of the condition on the average ratio of P t to P^ at the inner horizon required by conformal separability, 
equation (|160|) . to the constraint on the same ratio imposed by the condition that particles free-fall from outside the outer 
horizon, equation (|16ip . for the case of an uncharged black hole with angular momentum parameter a — 0.96M. . In the 
equatorial region cos 9 ~ 0, outgoing particles (positive P t and P,f,) satisfy the constraint but ingoing particles (negative Pt and 
P^) do not. Conversely in the polar regions | cos#| ~ 1, ingoing particles satisfy the constraint but outgoing particles do not. 



different dependences on latitude. While the inequality (I161|) may be satisfied at any latitude by either the ingoing 
or the outgoing stream, the inequality cannot be satisfied simultaneously by both ingoing and outgoing streams at 
that latitude. 



D. Dispersion of angular motions incident on the inner horizon 

The purely angular components of the collisionless energy-momentum are sub-sub-dominant. During inflation, the 
angular components are negligible because the incident accretion flow is by assumption tiny, and inflation does not 
amplify the angular components of the momenta of the streams. During collapse, the angular components of the 
momentum grow faster than the radial components, eventually becoming comparable to the radial components, at 
which point the solution fails. As argued in §VIII I[ if the Einstein equations are required to hold also for the sub- 
sub-dominant angular components of the collisionless energy-momentum, then conformal separability requires that 
the angular energy-momentum be isotropic (proportional to the unit 2x2 matrix) . If the angular energy-momentum 
is initially isotropic, then the behaviour of freely-falling streams ensures that it will remain so during inflation and 
collapse. 

An isotropic angular energy-momentum incident on the inner horizon can be contrived by allowing multiple ingoing 
and outgoing streams whose angular components of momentum satisfy condition (| 1 29|) in the mean, but are isotropic 
in their mean squares, which can be arranged. 

It is apparent that the required boundary conditions become more special and more contrived as the condition of 
conformal separability is tightened. To satisfy the radial components of the Einstein equations, only condition (|158[) 
on the incident densities N ± is required. To satisfy the sub-dominant radial- angular components of the Einstein 
equations, condition (|129j) on the angular components of the momenta is required. To satisfy the sub-sub-dominant 
angular components of the Einstein equations, the angular energy-momentum tensor must be isotropic. 

XI. CONCLUSIONS 

This paper has presented details of a new set of conformally stationary, axisymmetric, conformally separable 
solutions for the interior structure of a rotating uncharged black hole that accretes a collisionless fluid. The solutions 
are generalized to charged black holes in Paper 3 |2J. A concise derivation appears in Paper 1 [3j. 

The solutions confirm explicitly, for the special case of conformal separability, that inflation develops in rotating 
black holes as anticipated by [Hj]. 
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Conformal stationarity means that the black hole accretes steadily at an asymptotically tiny rate, in such a way 
that the spacetime expands conformally (self-similar ly), without change of shape. The solutions presented should 
be a good approximation, in the sense of perturbation theory, when the accretion rate is small, and should become 
increasingly accurate as the accretion rate tends to zero. 

The most important equations in this paper are equations (1881) for two of the diagonal Einstein components, which 
have negligible collisionless source of energy-momentum. These equations hold over the entire regime of interest, from 
electrovac through inflation and collapse. The equations can be separated, and their solution seamlessly yields both 
electrovac and inflationary solutions. The equations recover several known physical features of inflation: inflation 
occurs at the inner horizon but not at the outer horizon; inflation requires the simultaneous presence of both ingoing 
and outgoing streams; and the smaller the accretion rate, the more violently inflation exponentiates. 

A central feature of separable electrovac (A-Kerr-Newman) solutions is that as ingoing (p t < 0) and outgoing 
(pt > 0) collisionless streams approach the inner horizon, they concentrate into narrow, intense beams focussed along 
the ingoing and outgoing principal null directions, regardless of the initial angular motion of the streams. If there 
were no back-reaction on the electrovac geometry, then the streams would exceed the speed of light relative to each 
other and fall through distinct ingoing and outgoing inner horizons into two causally separated regions of spacetime. 
In reality, the energy and pressure of the counter-streaming ingoing and outgoing beams builds to the point that 
it becomes a significant source of gravity, however tiny the accretion rate. The gravity of the streaming energy- 
momentum acts so as to accelerate the ingoing and outgoing beams even faster through each other along the principal 
null directions. The result is the inflationary instability, in which the energy-momentum of the hyper-relativistically 
counter-streaming beams grows, along with the Weyl curvature and mass parameter, to exponentially huge values, 
while the radial horizon function decreases exponentially. 

At an exponentially tiny value, the horizon function goes through a minimum (in absolute value), whereupon the 
spacetime collapses, the conformal factor shrinking to an exponentially tiny scale. This is consistent with the argument 



of [11| that inflation leads to collapse, not a null singularity, if the black hole continues to accrete, as is ensured in 
the present solutions by the presumption of conformal time-translation invariance (self-similarity). 

During collapse, the angular motion of the infalling streams grows (the tetrad is chosen to align with the principal 
frame, so the angular directions are orthogonal to the principal null directions). The solutions presented here break 
down when the angular motion of the streams becomes comparable to the radial motion, but this occurs only when 
the spacetime has collapsed to an exponentially tiny scale. That the solutions fail when angular motions become large 
is consistent with the physical idea that conformal separability can persist only so long as the energy-momentum is 
focussed along the radial direction. The result is also consistent with the slowly rotating black hole solutions of [l7j . 

Until angular motions become important, the radial components of the collisionless energy-momentum dominate 
the angular components, satisfying the hierarchy of inequalities 

\T XX \ » \T xy \ » \T m \ , (162) 

where x signifies radial (x-t) directions, and y angular (y-(f>) directions. The radial Einstein equations are unaffected 
by the sub-dominant radial-angular Einstein equations, which in turn are unaffected by the sub-sub-dominant angular 
Einstein equations. 

Conformal separability imposes a corresponding hierarchy of boundary conditions on the collisionless ingoing and 
outgoing streams incident on the inner horizon. The indispensible boundary condition is set by requiring that the 
radial Einstein equations be satisfied. Conformal separability requires that the accretion flow be "monopole" in the 
sense that the densities of ingoing and outgoing streams must be independent of latitude, equation (|158l) . It makes 
physical sense that conformal separability would require this high degree of symmetry on the incident accretion flow. 
If the accretion rate were different at different latitudes, then the spacetime would collapse faster where the accretion 
rate is higher, destroying the symmetry. 

The evolution of the conformal factor and radial horizon function during inflation and collapse is determined by 
two small constant parameters u and v set by the incident densities of ingoing and outgoing streams, which are 
proportional to u =p v. Positivity of both densities, coupled with the requirement that the black hole expand as it 
accretes, imposes u > v > 0. The case v = corresponds to the stationary (homogeneous) approximation of [32T434] ] . 

If the sub-dominant radial-angular Einstein equations are required to be satisfied, then conformal separability 
requires that the angular components of the momenta of the incident ingoing and outgoing streams have Hamilton- 
Jacobi parameters satisfying equation (|129l) . 

A limitation of the required angular motions is that they cannot be accomplished by collisionless streams that fall 
freely from outside the outer horizon. Particles that fall from outside are necessarily ingoing at the outer horizon. 
This condition excludes half the phase space available to freely falling particles, and makes it impossible to fulfill the 
required angular conditions at the inner horizon. Thus, if the angular conditions are imposed, then the ingoing and 
outgoing streams must be regarded as being delivered ad hoc to just above the inner horizon. 
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If the sub-sub-dominant purely angular Einstein equations are required to be satisfied, then conformal separability 
requires that the mean squares of the angular components of the collisionless momenta be isotropic, which can be 
contrived. 
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Appendix A: Reduction of the line-element 

This Appendix shows how the assumptions of conformally stationarity, axisymmetry, and conformal separability 
lead to the line-element ([3]) and associated vierbein © adopted in the text. 

By conformal separability is meant conditions on the vierbein e m M and electromagnetic potential A m that follow 
from assuming that the action S governing the free motion of neutral or charged particles separates as a sum of terms 
each depending only on a single coordinate, equation (|28p. and that the left hand side of the resulting Hamilton- Jacobi 
equation, after multiplication by an arbitrary overall factor, separates "in the simplest possible way," assumption III 
of |3lj . as a sum of terms each depending only on either radius x or angle y. 

Conformal separability differs from separability in that the resulting Hamilton- Jacobi equation separates exactly 
only for massless particles, m = 0. Consequently the spacetime has a conformal Killing tensor. 

The assumptions of conformal stationarity and axisymmetry imply that the canonical momenta 7r t and ir^ are 
constants, equations (|22[) . The assumption of conformal separability implies that the canonical momenta n x — dS/dx 
and TTy — dS/dy are functions only of x and y respectively. 

Let e m ' 1 and A m denote the vierbein coefficients and tetrad-frame electromagnetic potential with an overall 
conformal factor p factored out: 

Cm'' = P&m^ ) A m — pA m . (Al) 

The conformal factor p here could be any arbitrary function of all the coordinates. In terms of the scaled vierbein 
e m M and electromagnetic potential A mi the Hamilton- Jacobi equation (|17[) is 

V" m (em M % - qAjj (e n »n u - gi n ) = -m 2 p 2 . (A2) 



Following assumption III of [311] ■ assume that the left hand side of equation (|A2I) separates in the simplest possible 
way as a sum of terms each of which depends only on x or only on y. Since 7r t and 7r^ are constants, while n x and tt v 
are respectively functions of x and y, the assumption implies that 

. . [ either L 11 for all u, and Am, are functions of x only, and e m y = , , . „, 

for each m, I m l ' i r r , , ( A3 ) 

^ or &7n f° r au Mi an d A m , are functions of y only, and e m — . 
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The case that matches the A-Kerr-Newman geometry, which provides the boundary conditions for the inflationary 
solutions considered in this paper and its companions, is 



the 



top 

bottom 



condition of (|A3j) holds for 



m = x and t 
m = y and <f> 



(A4) 



Thus conformal separability "in the simplest possible way" consistent with A-Kerr-Newman requires that 



v — £.v — p x — s x — 



= 



(A5) 



Given the conformal separability conditions (IA5|) . the vierbein coefficients it x and e$ v can be transformed to zero 
by a tetrad gauge transformation consisting of a Lorentz boost by velocity e t x /e x x between tetrad axes ~f x and f t , 
and a (commuting) spatial rotation by angle tan -1 (e^ / 'e y y ) between tetrad axes 7 y and 70. Thus without loss of 
generality, 



e t x = = 



(A6) 



The gauge conditions (|A6I) having been effected, the vierbein coefficients e x , e y *, e x ^, and e y ^ can be eliminated by 
coordinate gauge transformations t t' and cf> (f>' defined by 



dt = dt' + 



e * e 4 

'1 6x dx + ^ 
e y v 



dy , 



p p V 

ex 1 e,/ 



(A7) 



Equations (|A7[) are integrable because e x ^ and are respectively functions of x and y only. The transformations (|A7[) 
of i and </> are admissible because they preserve the Killing vectors <9/9i and d/dcj), 



a 




d 




d 


d 


dt 




~ dt' 


x,y,4> 


d(j) 


•A, ~ W 



(A8) 



x,t,y 



Thus without loss of generality 



6-r — 67, Gj 



e/ = 



Finally, coordinate transformations of the a; and ?/ coordinates 



y^y 



(A9) 
(A10) 



can be chosen such that e x x and e y y satisfy 



(All) 



The conformal separability conditions (|A3[) with (JA4J) , which imply conditions (|A5|) . coupled with the gauge 
conditions (| A6|) . (IA9|) . and (jAlip . lead to the line-element ([3]) and vierbein ((SJ) adopted in this paper. 



Appendix B: Tetrad-frame connections, Einstein and Weyl tensors 



This Appendix gives expressions for the tetrad-frame connections, and Einstein and Weyl tensors, in the case that 
the conformal separability conditions (|29[) hold, and the conformal factor p is any arbitrary function not only of x 
and y, but also of t and <j>. There arc 18 non- vanishing tetrad-frame connections T/c/m (counting — — ^ikm &s 
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one), comprising 8 distinct connections, 
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^txx 




= r^ t 
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= r t ,, t 


r (frxx 


— Lt(/>t 


= Ty^j, 


Fxt<j> 




= ^tjitX 




= r y f 


= 












^<j>y<j> 


^xy4> 




— ^ytpx 



d x In p , 
<9 t In p , 
<9 y In p , 
d$ In p , 

2p<r 2 



2pa 2 
^ In p 

<9 a In p - 
. 



p <9x 
a 2 dy/A^/cr 2 

p ay 



(Bla) 
(Bib) 
(Blc) 
(Bid) 

(Blc) 

(Blf) 

(Big) 

(Blh) 
(Bli) 



Since the tetrad-frame connections Tkim constitute a set of 4 bivectors, it is possible to combine the connections into 
complex combinations, but no additional insightful result emerges from those combinations. 
The 10 tetrad-frame components Gki of the Einstein tensor satisfy 



-D k dip + Gki 
P 



mi 



iR^-dmpd" 



= < 



4C s -2C v +rl t<p + 


1 yt<P 


XX 


-4C X + 2C V -J% 


4- r 2 


tt 


-AC y + 2C x + rl t<p 


+ r 2 

+ 1 yt4> 


yy 


-4C V + 2C X -I%^ 


+ r 2 
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— ^xt4^yt4> 




xy 


-2(W X - Wy) 




t<f> 







xt,x<f>,yt,y(f> 



(B2) 



where R is the Ricci scalar, and C x , C y , W x , and W y are related to the Weyl tensor as described immediately below. 
The Weyl tensor can be thought of as a matrix with bivector indices, and as such has a natural complex structure 
[39| . embodied in the complexified, self-dual Weyl tensor Ckimn defined in terms of the usual Weyl tensor Ckimn by 



a 



■£kf 



a 



pqrs 3 



(B3) 



where Sklmn is the totally antisymmetric tensor, normalized here to e klrnn = \klran] in an orthonormal tetrad frame. 
The distinct components of the complexified Weyl tensor form a 3 x 3 traceless, symmetric, complex matrix, whose 
5 distinct complex components form objects of spin 0, ±1, and ±2. The spin-0 part of the Weyl tensor is 



C C x t x t ^C X yxy — r^Cxtxcf) 



(B4) 



The real and imaginary parts of the spin-0 Weyl tensor are commonly called its polar (p) and axial (a) parts, the 
polar part remaining unchanged, while the axial part changes sign, under a flip 7^ — » — 7^ of the azimuthal tetrad 
axis: 



C = C {p) + iC {a) . 
In the present case, the polar spin-0 component of the Weyl tensor is 

C^ = C X +Cy , 

where 

Udx 

12 dy |_ dy 



P 2 C X 



p 2 Cy = 



cr 



2 d(A x / er 4 ) 



dx 
d(A v /a*) 



1 dui x 
a 1 dx 
1 dujy 
a 2 dy 



A, 



(B5) 
(B6) 

(B7a) 
(B7b) 



3G 



while the axial spin-0 component is given by 

few = - 



dujy d(A x /a i ) _ du x d(A y /a i ) 
dy dx dx dy 



The only other non- vanishing component of the Weyl tensor is the spin-1 component 

C xtX 4> = W X + W y , 

where 



2 W _ \/-&x^y d f I dui x 
4 dx \<t 2 dx 

„2 W _ ^J-^x^y d ( 1 dwy 
PWv ~ 4 dy U 2 dy 
In the spacetimes presented in this paper and its companions, the spin-1 component (IB9[) always vanishes. 



(B8) 

(B9) 

(BlOa) 
(BlOb) 



Appendix C: Electrovac solutions 



The standard electrovac solutions can be derived from the assumptions of strict stationarity, axisymmetry, and 
strict separability as follows. As shown by [3l|, the line-element takes the form ([3]) with a separable conformal factor 
p s , equation (|3"0)) . 

Given strict stationary (v = 0) and separability, the Einstein component G xy , which has zero electrovac source, is 



PsGxy 



;^-A x A, 



- d 2 Hp 2 J° 2 ) 

dxdy 



(CI) 



Homogeneous solution of this equation implies the form (|72p of the conformal factor p s . At this point the constants 
go and g\ can be adjusted arbitrarily without affecting either p x or p y : the overall normalization of go and g\ is 
cancelled by the normalizing factor of l/V/o.9i+/iffo m Px and p y , and the relative sizes of go and g\ can be changed 
by adjusting an arbitrary constant in the split between p x and p y . 

The Einstein component Gt$, which also has zero electrovac source, is 

p 2 G t </> - 



V-&xA y 




( Ps du x \ 




( pj dujyY 


2p 2 


dx 


\a 2 dx ) 


- — 

dy 


ya 2 dy J 



(C2) 



which, given the expression (f72j) for the conformal factor p s , reduces to 

1 y / -A x A y \du x dln(/ + f\ui x ) du/ y dln(/i + f u x ) 



P s G t 4> — - 



dx 



dx 



dy 



dy 



(C3) 



Homogeneous solution of this equation can be accomplished by separation of variables, setting each of the two terms 
inside square brackets, the first of which is a function only of x, while the second is a function only of y, to the same 
separation constant 2/2. The result is 



duj x 
dx 



2\ Ufo + fiWx) 



go + -rihgo + h)t^3 

Jo 



dui y 
dy 



= 2J(f 1+ f0LOy) 



91 + -rifogi + h)u v 
h 



(C4a) 
(C4b) 



for some constants go and g\ , which can be taken without loss of generality to equal those in the conformal factor ([72 
The Einstein components G xx + Gu and G yy — G^a,, which too have zero electrovac source, are 



pl (G xx +G tt ) = 



Ps (Gy 



Ga 



2A 
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a 2 




PsTT- 

dx 


2Ay 




9 ( 


a 2 


HS dy V 



•d(l/Ps) 



dx 

2 9(1/P S ) 



1 



dujo 



dy 



4cr 2 \ dy 

1 / dux 
4cr 2 V dx 



(C5a) 
(C5b) 
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which with the conformal factor p a given by equation (|72[) and dui x /dx and duj y /dy given by equations (|C4[) reduce 
to 



P 2 S (G xx + G tt ) = -2A 

Ps ( G yy ~ G 00) = 2 \ 
These vanish provided that the constant ji satisfies 



{h + /offi + /igo)(/o + fiu x ) 2 
/0M 4 

(/ 2 + /o<?i + /i<?o)(/i +/ow s ) 2 



fofi<? 4 



(C6a) 
(C6b) 

(C7) 

Inserting this value into equations (|C4|I yields equations ([73]) . 

The four off-diagonal Einstein components G^t, G y t, G x ^, and Gj,^ vanish identically. 

Solution of the Einstein equations for the remaining two Einstein components G xx — Gu and G yy + G^ was already 
discussed in ^VIII Dl Inserting the homogeneous solutions (|90c[) and (|90dl) into the defining equations (|87cl) and (|87dp 



/a = -(/offi + /iffo) • 



for Y x and Yy, and setting U x = U y = (which eliminates inflation), yields differential equations for the radial and 
angular horizon functions A x and A y , solution of which, with appropriate boundary conditions, recovers the Kerr 
solution. 

Solutions including the energy-momentum of a static electromagnetic field fall out with little extra work. With 
appropriate boundary conditions, this is the Kerr-Newman solution. Electrovac solutions have Gtt — —G xx and 
G^ — G yy . Electrovac solutions with G yy = G xx , as is true for a static radial electromagnetic field, are found by 
taking the difference of equations (|88| and separating that difference in the pattern of equation (|89l) . The solution is 
a sum of a homogeneous solution and a particular solution 



Y, 



HQ 2 . + QlKfo + fiu*) 

dui x /dx 



Y u = 



(C8) 



Inserting equations (|C8[) into the Einstein expressions (|88|) yields Einstein components that have precisely the form 
Gmn = [(Ql + Ql)/pf\ diag(l, -1, 1, 1) of the tetrad-frame energy-momentum tensor of a static radial electromagnetic 
field. 

Similarly, solutions including vacuum energy, which has G yy = —G xx , can be found by separating the sum of 
equations (|88p in the pattern of equation 



Yr 



A particular solution is 
2A 



fldu x /dx 



Y k 



2Aw 2 



ffduj x /dx 



(C9) 



Inserting equations (|C9|) into the Einstein expressions (|88|) yields Einstein components that have precisely the form 
of a cosmological constant, G mn = —Arj mn . 

Solving equations (|87c|) and (|87d[) with Y x and Y y given by a sum of the homogeneous, equations (|90cj) and (|90d[) . 
electromagnetic, equation (|C8I) . and vacuum, equation (IC9[) . contributions, yields the general electrovac solution for 
the radial and angular horizon functions A x and A y , 



= (fa + fiw x ) 



&y = (fl + foUy) 



(k + kiu w ) - 



(fci + koLo y ) 



2M.y/(f + f 1 Lj x )(g -9iu x ) , (Q5 + Q2)(/o + /iWx) 



(/o.9i + /1.90) 3 / 2 

2A/'.- v /(/i + fou y )(gi - 90Uy) 
(fo9i + /150) 3 / 2 



/0.91 + /150 
3/i (fo9i + /i.9o) 2 



A(.9o - 9iu x ) 
3/i(/offi + /iffo) 2 



(ClOa) 
(ClOb) 



where fc and fci are arbitrarily adjustable constants arising from the freedom of choice in the constants h and hi of 
the homogeneous solution. The constant M, in the expression (IClOap for A x is the black hole's mass. The constant 
Af, in the expression (jClObp for A y is the NUT parameter 0, |40I442| |. which is to the mass M, as magnetic charge Q, 
is to electric charge Q». 

The electrovac Weyl tensor (|B5| has only a spin-0 component, and is 



C 



(Px - iPyf 



M. + iAf. 



Ql + Ql 

Px + iPy 



(Cll) 



with p x and p y given by equations (|72bl) . 
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Appendix D: Integrals along the path of a particle 



This Appendix derives conditions under which integrals along the path of a particle can be deemed small, in the 
conformally stationary limit of small accretion rate. The results confirm that the tetrad-frame momentum p^. and the 
density N predicted by the Hamilton- Jacobi equations for a massive particle are accurate, by demonstrating that the 
integrals of equations (|39|) and (j63j) are adequately small. 

Consider an integral of the form 

" \A x \ a (p/p s fdX . (Dl) 

Along the path of a particle, dX — (p 2 /\P x \) dx, and p/p s — e vt ~^ — [(u=fv)/(U x 1= d)] 1 ' 2 from equations (|M)l and 
(|114p . The integration interval dx may be converted to either dA x or dU x using either of equations (|93l) . The 
integral (|Dlj) thus becomes 

|Ax| 1+Q ( u T v \ 1+P/2 d\n\A x \ f {A^ ( u T v \ 1+p/2 dU x 



P*\ \U x TvJ A' X -3U X J \P X \ \U x tvJ 2(U 2 -v 2 ) 



The horizon function A x is given as a function of U x by equation (|99|) . 

Inflation ignites when |A X | ~ v, and the integral (|D1[) may be counted from this point. At the beginning of inflation, 
before lA^I has been driven to an exponentially small value, U x changes little from its initial value of it, and then the 
middle expression of equations (|D2[) yields 

IA l 1+a 

This integral is small in the conformally stationary limit v — > provided that a > — 1, but as to whether an integral 
in an actual situation can be judged small will depend on additional factors of v that the integral may be multiplied 
by 

By the time that U x has reached of order unity times its initial value u, the horizon function \A X \ has become 
exponentially small. In this regime further contribution to the integral is exponentially small provided that 

a > -1 . (D4) 

During collapse, both lA^I and U x increase exponentially (|A X | from an exponentially small starting point). As 
long as |A X | < 1, further contribution to the integral remains exponentially small provided that 

/3 > -4 . (D5) 



At the end of collapse, there is a regime where \A X \ > 1 and \P X \ ~ \/—A x but |A X | <C U x , while U x is exponentially 
huge. Here further contribution to the integral remains exponentially small provided that 

either /3 > -3 + 2a or f3 > -4 . (D6) 

In the case of equation (|31?|) for dpk/dX there are two terms. The first has a — and (3 = —3. The second has a = 
and P — — 1, except that the x component has a = —1/2 and f3 = — 1 when lA^I < 1. Both terms have pre- factors 
proportional to v during early inflation (including the term proportional to d{p 2 — p 2 ) dy). For both terms, the main 
contribution to the integral is from early inflation, equation (|D3|) . which, multiplied by the v pre- factor yields ~ v 2 
for the first term, and ~ v 3 ^ 2 for the second. The ratio of the integral of dpk/dX to pk itself is, since p x oc 1/ y / jA7f 
when |A X | < 1, 

Apfe _ J dp k /dXdX 2 



Pk Pk 



v 2 . (D7) 

Both terms satisfy the conditions (|D4|) . (|D5j) . and ()D6|) so the contributions to the integrals after early inflation are 
exponentially small. 

In the case of equation (j63]) for Dkn k , the exponents are a = and /3 = — 2, and the pre- factor is proportional to 
v during early inflation. Again, the main contribution to the integral is during early inflation, giving 

Ainjv _ [dpgdXdX ^ y2 
Pk 

Again the conditions (|D4I) . (|D5I) . and (|D6|) are satisfied, so the contribution to the integral after early inflation is 
exponentially small. 



